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Abstract—Carbon dioxide was absorbed into water containing surface-active agenis and 
running in a thin film over a sphere. 

The results are compared with theoretical calculations based on the assumption of equilibrium 
at the interface, absorption being controlled by diffusion. The comparison indicates that at a 
certain concentration of the commercial agents there is a resistance to absorption which falls 
to zero at very low and at very high concentrations. With a chemically pure agent, there is no 
resistance to absorption, other than that due to diffusion, so that theory and experiment agree. 
With a solution containing the same pure agent plus a small amount of impurity, the behaviour is 


similar to that of the commercial agents. 


These results indicate that surface-active agents can reduce the absorption up to 25% in 
a system in which they do not influence the hydrodynamics. 


Résumé 
s'écoulant en couche mince sur une sphére. 


De l’'anhydride carbonique est absorbé par de l'eau contenant des agents tensio-actifs, 


Les auteurs comparent les résultats obtenus avec les calculs théoriques basés sur I'hypothése 


d'équilibre 4 linterface, la diffusion contrélant l'absorption. 


La comparaison montre que pour 


une certaine concentration des agents commerciaux, il se produit une résistance a l'absorption 
qui tombe a zéro pour des concentrations trés faibles ou trés élevées. 

Avec un agent chimiquement pur, il n'y a, comme résistance a l'absorption, que celle inhérente 
4 la diffusion, ce qui explique l'accord entre la théorie et l'expérience. Une solution contenant lé 
méme agent pur, avec une petite quantité d'impuretés, se comporte comme les agents commer- 


ciaux. 


Ces résultats montrent que les agents tensio-actifs peuvent réduire l’'absorption jusqu ’& 25% 
dans un systéme oi ils n’influencent pas [hydrodynamique. 


INTRODUCTION 

Tue effect of surface active agents on such 
properties of an interface as the surface tension 
is well known, but their influence on the rate 
of mass transfer across an interface is uncertain. 
This uncertainty is due to the fact that mass- 
transfer experiments are usually carried out on 
moving systems, and the addition of surface- 
active agents may alter the hydrodynamics, 
an effect which is difficult to separate from a 
surface resistance. 

In liquid-liquid systems, two methods have 


been used. In the first of these, the effect on the 
rate of mass transfer to or from a drop, which is 
rising or falling through a stagnant liquid, is 
measured. GarNeR and his co-workers [1, 2] 
have used this method and concluded that surface 
active agents present a resistance to mass transfer. 
LinDLAND and TERJESEN [3] came to a similar 
conclusion and developed a theory for the resis- 
tance on the basis that it is proportional to the 
interfacial concentration of surface-active mole- 
cules. The second method is to measure the 
transfer rate across an interface between stirred 
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liquids. Hutcninson [4] concluded from this 
type of experiment that sodium cetyl! sulphate 
caused an interfacial resistance, the magnitude 
of which was dependent on the diffusing molecules. 
Hoim and Tersesen [5] found a minimum 
in the transfer rate and concluded that the rise 
after the minimum was due to an increased 
interfacial area. An excellent review of earlier 
work on these two methods and on the effect 
of surface films on evaporation is given in reference 
{1}. 

The data in the literature on the effect of surface 
active agents on the rate of absorption of gases 
by liquids is not so extensive. TeRNOVSKAYA 
and Be.opo.skii [6] studied the effect of three 
sulphonate surface-active agents on the rate 
of absorption of sulphur dioxide in water flowing 
down a wetted-wall column. With increasing 
concentration of one of the agents, whilst the 
surface tension decreased, the HTU 
at first, and then fell off after reaching a maximum 
at a certain concentration. The effect of the 
second agent was similar but less pronounced, 
while the third caused an HTU 
to a constant value. These effects were said to 
be due to the adsorption layer at the interface 
modifying the hydrodynamics [7, 8, 9]. EmmMert 
and Picrorp [10] compared the theoretical and 
experimental rates of absorption and desorption, 
for oxygen and carbon dioxide, in water which 
a 


increased 


increase in up 


contained sulphonate surface-active agent, 


in laminar motion down 
wetted-wall column. They concluded that there 
interfacial due to lack of 
equilibrium at the interface, even with pure water, 
and that the sole effect of the agent is to eliminate 
rippling. Subsequent work done by the present 
authors [11, 12] suggests that equilibrium does 
exist at the surface of pure water, and that 
the discrepancy between the theoretical and 
observed results obtained by Em™Mert and 
Picrorp is due to an interfacial resistance caused 


by the surface-active agent. 


and which flowed a 


is an resistance, 


EXPERIMENTAL 


In the present investigation, a modified form 
of wetted wall column has been This 
column, which utilizes water flowing over a 


used. 


50 


sphere instead of down a cylinder, has been 
described in detail elsewhere [11]. 

The absorption gas used was carbon dioxide, 
produced at a purity greater than 99.8%, from the 
solid. Air-free distilled water, with varying 
concentrations of surface-active agents, was 
used as the absorbing liquid and its temperature 
was maintained at 25-04 1°C. The whole 
apparatus was enclosed in a constant-temperature 


room which was kept at 25.0 + 05°C. During 


an experiment the carbon dioxide was drawn 


in from an atmospheric balloon through a rota- 
meter which measured the rate of absorption. 
The column was purged at frequent intervals 
to prevent the buildup of impurities. 

Four commercial surface-active agents were 
used, together with a specially prepared sample 


of one of them. These were : 


Teepol 


21%, active secondary alkyl sulphate with 


An aqueous solution containing 


average molecular weight of 305, and a 
specific gravity of 1.075 at 25°C. 
R 


taining approximately 33° 


Petrowet An aqueous solution con- 


» of saturated 
hydrocarbon sulphate with an average 
molecular weight of 280 and a specific 
gravity of 1-067 at 21-5°C, 
Lissapol N — An 


taining 265°, of ethylene oxide con- 


aqueous solution con- 
densates of an average molecular weight 
of 638, and a specific gravity of 1-024 
at 21-5°C [13]. 


(4) Sodium dodecyl A solid 


commercial product of BP purity. 


sulphate 


A specially pure sample of sodium dodecyl 
sulphate was prepared from homologue-free 
dodecyl alcohol. The alcohol was sulphated 
with chlorosulphonic acid in the presence of 
acetic acid at O°C; after neutralizing with 
sodium carbonate, the sulphate was extracted 
with butanol, and after evaporation, crystallized. 
The crystals were thoroughly dried and were 
extracted with ether in a Soxlet for 16 hours. 
The product was then crystallized from ethanol 
and carefully dried. 

For the liquid agents, solutions were made 





VOL. 
6 


956/57 


The effect of surface active agents on the rate of absorption of carbon dioxide by water 


up to a known volume percentage and all con- 
centrations referred to in this paper are with 
reference to the material as supplied and not 
to the active component. For the solid agents, 
solutions were made up to a known weight 
percentage of the material as supplied or made. 

For ease of comparison, all experimental 
absorption rates were corrected to a partial 
pressure of carbon dioxide of 760 mm and to a 
water temperature of 25°C. The procedure 
for doing this is given in reference [11] and the 
maximum total correction needed in this work 
was under 8%. 

To characterize the solutions, measurements 
were made of the surface tension and equivalent 
conductivity over the concentration range 
required. The surface tensions were determined 
using the drop weight method, with the corrections 
of Harkins and Brown [14]. The electrical 
conductivity was measured using a Mullard 
conductance bridge and a conductivity cell. 
The conductivity of the agent was found by 
subtracting the conductivity of the distilled 
water (about 38 x 10% ohm ™'cm™"') from the 
conductivity of the solution. With solutions 
of Teepol and Petrowet, which are mixtures of 


several substances of varying equivalent weights, 


the equivalent conductivity based on an average 
equivalent weight has been used. The presence 
of inorganic salts and other compounds added 
to enhance the commercial value of these agents 
renders the absolute value of these “ equivalent 
conductivities ” doubtful. 


RESULTS 

The effect of various concentrations of the 
surface active agents on the rate of absorption 
of carbon dioxide by water flowing over the sphere 
was determined. Flow rates between 1 and 
7-5 cc/sec of water were investigated and con- 
centrations up to a maximum of 5%, by volume. 
As shown by Figs. 1-4, the results for the 
commercial products were similar to one another 
in that, at a given liquid rate, the absorption 
fell sharply with small concentrations to a 
minimum, and eventually at much larger con- 
centrations returned nearly to the theoretical 
rate. The theoretical rate was calculated as 


CO, absorbed G, g/sec x 10* 
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Fic. 1. Absorption rates with Teepol solutions. 
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Fic. 2. Absorption rates with Petrowet solutions. 
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Fic. 4. Absorption rates with commercial sodium dodecyl] 
sulphate solutions. 


in reference [11]. With very low concentrations, 
and at liquid flow rates greater than 3 cc/sec, 
the absorption rates were greater than the 
theoretical due to rippling in the liquid film. 
The rippling was eliminated when the wetting 
agent exceeded a certain concentration, this 
being less than the concentration to give the 
minimum absorption rate. 

Fig. 5 gives results for the pure sodium dodecyl 
sulphate. A solution containing 0-1 wt.°, was 
used over the whole range of liquid flow rates, 
and at a flow rate of about 2 cc, sec, concentra- 
tions of 0-08, 0-2 and 0.3 wt.°®,, were investigated. 
The results, which are compared with the theore- 


tical absorption rate, show that pure sodium 


dodecyl sulphate has a negligible effect on the 
absorption. 

Fig. 6 shows the variation of the absorption 
rate, with concentration, for a solution of pure 
sodium dodecy! sulphate plus 2°,, dodecanol (based 
on sulphate), at a liquid flow rate of 2 cc /sec. 

Figs. 7 and 8 show the variation of the surface 
tension and equivalent conductivity with volume 
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Absorption rates with pure sodium dodecyl 
sulphate solutions. 
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Weight percentage 
Fic. 6. Variation of absorption rate with concentration of 


pure sodium dodecyl sulphate + 2°, dodecanol solutions. 
The horizontal dotted line represents the theoretical rate. 


concentration of Teepol and Petrowet R, together 
with the variation of the rate of gas absorption 
at liquid flow rates of 2 and 5cc/sec. Fig. 9, 
for Lissapol N, is similar except that, being 
non-ionic, no conductivity data are recorded. 
Fig. 10, for commercial sodium dodecyl! sulphate, 
is also similar, except that it is based on a weight 
instead of a volume percentage. These figures show 
that there is a rough connexion between the change 
in gas absorption and the changes in surface 
tension and conductivity. At large concentrations 
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Fic. 7. Variation of absorption rate, surface tension and 
equivalent conductivity with concentration of Teepol. 
The horizontal dotted lines represent the theoretical 
absorption rate. 
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Volume percentage of petrowet 
Fic. 8. Variation of absorption rate, surface tension and 
equivalent conductivity with concentration of Petrowet. 


The horizontal dotted lines represent the theoretical 
absorption rate. 


the absorption rate returns very nearly to the 
theoretical value. The variation of surface 
tension and equivalent conductivity with concen- 
tration for the pure sodium dodecyl sulphate 
was in good agreement with the values in the 
literature [15, 16]. All the surface tension 
and conductivity experiments were carried out 
near room temperature and the exact tempera- 
tures are recorded on the individual curves. 
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Fic. 9. Variation of absorption rate and surface tension 
with concentration of Lissapol. The horizontal dotted 
lines represent the theoretical absorption rate. 
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Fic. 10. 
and conductivity 
sodium dodecyl sulphate. The 

represent the theoretical absorption rate 


Variation of absorption rate, surface tension 
with 
horizontal 


equivalent concentration of 


dotted lines 


DiscussIoNn 


The results shown in Figs. 7 to 10, using the 
commercial surface active agents, show that the 
region of depressed absorption rate is about the 
same as the region in which the surface tension 
and conductivity are changing. These changes 
are generally attributed to the formation of 
micelles of the surface active molecules or ions 
in the bulk of the solution. As the concentration 
this the 


resistance to gas absorption gradually decreases 


is increased beyond critical value, 


until at relatively high values it becomes 
negligible. 

The suggested explanation of this effect is 
the that MILeEs 


Suep.ovsky [17] for a related effect. They found 


same as proposed by and 


that small quantities of dodecanol added to 
solutions of very pure sodium dodecy! sulphate 
caused a minimum in the graph of surface tension 
This was attributed to 


versus concentration. 


the sparingly soluble dodecanol being preferen- 
tially adsorbed at the surface at low concentra 
tions, but being desorbed back into the bulk 
of the solution after micelle formation had begun. 


The the commercial 
might well cause a similar effect on the rate of gas 


absorption. Support for this idea was obtained 


impurities in products 


from the results of experiments using very pure 
sodium dodecyl suphate, shown in Fig. 5. The 
The latter showed no resistance to gas absorption 
at any concentration up to 0.3 wt.°,. However, 
as shown in Fig. 6, the addition of 2°, of dodecanol 
(based on the sulphate) caused a high resistance 
to develop at concentrations less than 0-1 wt.°,. 

The reason for small quantities of impurities 
causing a resistance to gas absorption is very 
difficult to The that the 
impurities might cause a high viscosity in the 


find, possibility 
surface film, resulting either in a low diffusion 
coefficient or in a hydrodynamic effect, was 
investigated. Qualitative experiments were made 
on the movement of small particles on the surface 
of the solutions used, but only the impure sodium 
dodecy! sulphate and the pure material containing 
dodecanol showed a marked surface viscosity. 
In both cases, viscosity appeared to reach a 
maximum at a concentration corresponding to the 
maximum resistance to gas absorption and to 
This 


quantitatively by 


decrease with concentration. 


effect 


Durnam and Camp [18], who observed the same 


increasing 


has been measured 
general effect for solutions of pure sodium dodecy! 


sulphate containing dodecanol. In another 
system, Buakey and Lawrence [19] found that 
small quantities of a soap produced abnormally 
high surface viscosities in an alcohol film, but that 
the effect 


concentrations 


with increasing 


The 


showed no 


decreased to 
the 
however. 


zeTo 


above critical. other 


commercial products, 
appreciable surface viscosity, and thus it seems 
that the 
the surface 
transfer with impure sodium dodecyl sulphrte 


rough correlation observed between 


viscosity and resistance to mass 


is coincidental. However, more quantitative 
work on surface viscosities is required before 
definite conclusions may be drawn. 

Another possible cause of the resistance is an 
interaction between the adsorbed film and the 
diffusing This 
likely to be simply mechanical obstruction, owing 
to the of the 


The present experiments indicate that either the 


molecules. interaction is not 


extreme thinness surface film. 


impurities alone or some combination of the 
impurities and the surface active agent constitute 


the resistance to mass transfer. That the latter 
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is more likely is suggested by the results with 
pure sodium dodecyl sulphate, though the work 
of LInDLAND TeRJESEN [4] 
this idea. With increasing concentration of a pure 


and contradicts 
surface-active agent, they found a reduction to a 
constant value in the overall extraction coefficient 
of iodine from water by drops of carbon tetra- 
chloride. This attributed to a 
resistance, but it could equally be interpreted as a 
hydrodynamic effect in the 
circulation of the drop is suppressed by the 
surface active agent. This idea is supported by 
the sharp reduction of droplet velocities with 
small additions of the agent. Garner and Hae 
[1], using Teepol, measured the overall coefficient 
for transfer from drops and found that the 
minimum 
with 

GARNER 


was surface 


which internal 


coefficient decreased to a and sub- 


increased steadily increasing 
the 


SKELLAND [2], using more specific agents, found a 


sequently 
concentration of agent. and 
similar effect, and the rise in the coefficient after 
the minimum attributed to 
deformation and oscillation of the drops. 


was increasing 
How- 
ever, in the light of the present results, it seems 
probable that the rise is due to desorption from 
the interface of the impurities in the surface- 


active agents used. Thus the evidence on whether 


the impurities alone or some combination of the 


impurities and the agent produce the barrier 
to mass transfer is not conclusive, but the present 
work suggests it is the latter effect. 

The 
Bevoroiskit [6] on absorption in a wetted-wall 
column is more readily explicable. The initial 
increase in HTU on of all three 
agents is almost certainly due to the elimination 


data obtained by TerNovskaya and 


the addition 
of rippling in the liquid film. The subsequent 
decrease in HTU observed with two of the agents 
is probably due to the desorption from the surface 
of impurities. EmMMERT PiGrorD’s 
[10] on a wetted-wall column is similar. Their 


and data 
results showed a surface resistance close to the 
value determined in the present work for the 
same agent. Their suggestion that the decrease 
in HTU after the the 
reappearance of rippling in the liquid film was 
qualitatively investigated, using a number of 
agents in a long wetted-wall column. In no 


maximum was due to 


case was it possible to detect any reappearance 
of ripples even at concentrations up to 5 vol.%. 

SHerwoop and Ho tioway [20] found that 
surface active agents had a similar effect on the 
rate of desorption of carbon dioxide in a packed 
tower. The explanation given for the data on 
wetted-wall columns is directly applicable, but 
uncertainty about the effect of the agent on the 
degree of wetting gf the packing renders this 
further evidence less definite. 


CONCLUSIONS 


1. Impurities in commercial surface-active 
agents can cause a reduction up to 25% in the 
This effect reaches a 


for 


rate of gas absorption. 
maximum critical 
micelle formation and, with increasing concen- 
tration, finally decreases to negligible proportions, 
The magnitude of the effect is specific to the 


near the concentration 


material used. 


2. Solutions of very pure sodium dodecyl 
sulphate show no resistance to gas absorption, 
but the addition of dodecanol causes a resistance to 
absorption. 


3. The use of commercial surface-active agents 
to eliminate rippling may introduce a resistance 
to mass transfer. However, with very pure 
concentration of 


very 


materials or with a sufficient 


commercial material, this resistance is 


small. 


4. The modified wetted-wall column used in 
the present work is very suitable for investigating 
the effect of 
transfer, because the hydrodynamic conditions 


surface-active agents on mass 
are well defined at the surface and it is possible 


to get good agreement between theory and 
experiment without adding any agent at all. 
However, relatively large quantities of material 


are required. 
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Critical and terminal velocities in fluidization 


P. H. Pincupeck and F. Popper 


Coal Tar Research Association, Gomersal, nr. Leeds 


Abstract 


Critical fluidization and terminal velocities of spherical and nonspherical particles were 


determined in nitrogen. The ratios of these velocities, or of their particle Reynolds numbers, are 


shown to be function of the dimensionless group 


yw? 


it wre 


The ractical implications of the results are discussed. 


Résumé—Les auteurs ont déterminé dans l'azote les vitesses critiques de fluidisation et les 


vitesses limites pour des particules sphériques ou non sphériques. 
Les rapports de ces vitesses ou de leur nombre de ReyYNoLps correspondant 4a la particule 


apparaissaient étre une fonction du groupe sans dimension ci-dessus. 
Examen des conséquences pratiques de ces résultats. 


1 IntrropuctTion 


Tue application of the technique of fluidization 
in the chemical industry is growing rapidly 
and is far outstripping the fundamental knowledge 
of the physical processes involved. 

One aspect which has received much attention 
is that of the critical velocity at which fluidization 
starts and outstanding in this respect are the 
pioneering papers by Leva and collaborators [1], 
Witueim and Kwavk [2] and by van HeerpDEeN 
[3]. Less attention has been given, as far as 
fluidization is concerned, to the velocity at which 
particles are removed overhead from the bed, 
though this aspect has been studied in connection 
with such subjects as transport of solids. The 
work carried out in relation to fluidization has 
recently been reviewed [4] and terminal velocities 
in the fluidization of particles with water have 
received some attention by WiLnetm and 
Kwavuk [2]. Van Heerpen et al [5], from their 
work on heat transfer in fluid beds, obtained 
approximate ratios of terminal to critical mass 
velocities for fully turbulent and for stream- 
lined conditions. Leva [6] has studied the rate 
of removal of fine particles from a fluidized bed 
of fine and coarse components and has expressed 
his results in terms of a first-order rate equation. 
A useful summary unrelated to the subject of 
fluidization is also given by DaLLAvALLe [7]. 


A general relation of initial fluidization velocity 
and terminal velocity would be of interest to the 
designer of fluidized solids reactors, where it is 
desired to avoid overhead loss of particles above 
any particular diameter from the reactor. The 
present study is the outcome of experiments 
in fluidized catalysis, where it was necessary to 
prevent overhead loss of particles without use of 
auxiliary equipment but it is believed that the 
relations derived below should also assist in the 
design of filters and cyclones in conventional 
fluidized reactors. 


2. THEORETICAL 


The free fall of small particles above about one 
micron diameter through a fluid is expressed by 
Srokes’ Law, 

Re, — PolPs — Po) gd° (1) 
18 p? 

For larger particles this equation is not valid. 
A relation for all sizes of particles has been 
derived by Rusey [8] from theoretical considera- 
tions. The total force holding a particle in sus- 
pension at a constant height in a fluid current 
may be considered as the sum of viscous resistance 
and fluid impact 


nd® 


6 


? ; 
(Pp — Po) 8 = Bad pry + Pe Ur (2) 
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Converting this equation by substitution of 


TF for vr and solving the quadratic equation 
Po 
we have 


‘ 3 
Rey = —6 + | + 2 (Pe am Po) Po 4 (3) 
3 p? 


A plot of Re, vs d from the relation approaches 
the Stokes curve from equation (1) asymptotically 
as particle diameter decreases. 

A measure of the interval in which particles 
are fluidized, but not removed from a vertical 
vessel is the ratio of Rey Reg. 
decide on the equation for the critical particle 
Reynolds number (Rey, at which fluidization 
starts) to be used. Both Leva’s [1] and van 
Heexpen’s [3] based the 
the [9] but differ 
numerically. It is shown later that the results 
of the present work are fitted much more closely 
by the van Heerpen equation [3] which for 


It remains to 


equations 
CarMAN-Kozeny 


are on 


relation 


spherical particles reduces to 


0-00123 x 0-594 po (p, — po) g 4,” 


fo ; 


im 
if fo is not neglected. 

Dividing equation (3) by equation (4) and 
replacing d by d, (as defined by van Heerpen) 
for the generalized case for particles of any shape 

3 
Po (Ps af Po) d, & 
2 


~ : 
+ | 6-744 » 10"(- —— 
J Po (Ps — Po) i) 


p? 
Po) a? g 


This ratio is obviously a function of the dimen- 


- 82122 « 10° 


10° 
Po (Pp, 


+ 12489 


sionless group 
2 
pe 


Po (p, Po) d, & 


may call the fluidization number, 


In this form the equation is applicable to 


which 
Fil. 


all fluids and spherical solids and in the case of 


we 


gases or vapours where py may be neglected 


in compariscn with p,, the group may be reduced 


we 


PoPede B 1 
generalized form, is equivalent to —~ — 
2K 4 p 
by Witnerm and Kwavuk [2] to express 
generalized pressure drop-velocity relationships 
for all particles and fluids. 
The accuracy of relation (5) 
experimentally. 


to This dimensionless group, in its 


as used 


was tested 


3. EXPERIMENTAL 


Experiments were carried out in a 2-inch diameter 
glass tube with nitrogen at room temperature. For the 
largest particles a 3-inch tube was used to determine 
Rey and a l-inch tube for Rey. 

Reg was measured by finding the velocity at which the 
pressure drop across the bed equalled the static pressure. 
The terminal velocity for spheres was determined by 
gradually increasing the nitrogen flow till particles began 
to be elutriated at the top of the tube. The particles were 
collected and the largest diameter of the particles carried 
over was determined microscopically. 

With nonspherical partiles nitrogen was passed 
through the bed at such a rate as to remove the finest 
material and the bed maintained for some time in this 
condition. When apparently no more particles were 
elutriated at this rate the nitrogen velocity was increased 
slightly and the particles now carried over were collected 
and d,, the effective diameter, determined using vAN 
Hererpen’s technique [3]. 

The densities of the particles in mercury p,, were 
determined in a modified density bottle at room tempera- 
ture. Samples of the materials were sieved through 
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Critical and terminal velocities in fluidization 


Table 1. Fluidization in nitrogen at 18°C py = 1-18kg/m® pp = 1-78 x 10° kg/m sec. 





Material Shape Density Average Effective Go 
kg /m® diameter diameter kg/m? sec 
d, d, 
x 108 x 108 


Spherical 0-0126 0-079 
“Ballotini” | 0-0303 0-290 
0-055 0-705 
0-139 2-850 
0-417 
1-060 


Microspherical |§ Mainly ‘ 0-038 
catalyst spherical 33: 0-083 1-610 
0-168 4-660 

Microspherical Mainly 5 0-0013 0-0041 

VOL, alumina spherical 0-0040 0-0222 
6 0-022 0-286 
956/57 3 — —~ = — 
Ground Non- 5i 0-031 0-460 

silica spherical 38: 0-065 1-440 

0-135 4-300 








standard sieves selected so that the bulk had a spread of 
about 10% from the mean diameter. In the case of 
spherical particles the average diameter was determined 
by microscopic measurement of 50 to 100 particles. 

The average effective diameter d, of nonspherical 
particles could be obtained from the straight line relation- 
ship of d, and d,, the mean sieve diameter, which was 
found in a subsidiary experiment to exist for these particles 
(see Fig. 1). 
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4. RESULTS 








The properties of the solid particles used in the 
experiments are given in Table 1 together with 
the experimental results for the critical Reynolds 
number. In Fig. 2 the van HeEERDEN and Leva 
equations are superimposed on a plot of the 
experimental results for glass spheres and it is 
clear that the former results in a much better fit. 
Above a diameter of 5-5 x 10*m, Re, is no 
longer a function of d® and van HEERDEN’s 
equation is no longer obeyed. 

The terminal velocities for the particles are 
listed in Table 2. As the terminal velocity 
Particle diameter was always that of the smallest particle in the 

Fic. 2. particular size range the d, values in Tables 1 
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Critical and terminal velocities in fluidization 


and 2 do not correspond. The values for Re, 
given in Table 2 were obtained from plots of 


Re, vs. d, for each type of material. The ratio 


Re 
T found from the interpolated values of Ry 
9 
and the experimental values of Re, is also given 
in Table 2. The experimental values of Rey 
for glass spheres are plotted in Fig. 3 with 
Rusey’s equation. 
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104 10-3 


Particle Diameter 


10-2 


Experimental results (Glass spheres) 
Theoretical curve (Rubey) 


5. Discussion 

The equation for the terminal velocity used 
in this study (eq. 3) represents the experimental 
results reasonably well. The scatter about the 
line is probably largely due to experimental 
error, but some real deviation from RuBery’s 
equation seems to exist for the larger particles. 
For the present purposes, however, RuBry’s 
equation is a good representation of the results. 


The experimental critical Reynolds numbers 
with values below 10 are all slightly lower than 
the values calculated from the van HEERDEN 
equation, but a sufficiently good fit for the 
vAN HEERDEN equation is to be accepted. 

For the larger particles (see Fig. 2), however, 
According 
to vAN HeEeRpeEN the equation is only valid for 
up to Reynolds 
present 


the equation is obviously not valid. 


the region of viscous flow, ie. 
about 10. The 
indicate that the limiting value is correct. 


numbers of results 
Whereas in the region of viscous flow Re, 
the 
apparently changes gradually to 


is proportional to d} proportionality 
the 


This change in propor- 


_ as 
particle size increases. 
tionality occurs also in the terminal Reynolds 
number (dependence of d3 here also ceases 
10) and is formalised in Rusey’s 
A modification of the van HEERDEN 


around Rey; = 
equation, 
equation along similar lines should provide a 
relation valid for all conditions ; however, such 
an expression would prove very unwieldy and 
would, in any case, be of only limited usefulness. 

The ratios vr from the experimental results 

€o 

obtained in this study are plotted in Fig. 4 
together with the line calculated from equation 
(5). 


obtained 


Also included in this figure are values 
from the WILHELM 
Kwak [2] for fluidization of particles in water. 


data of and 


The fit for the spherical particles is very good. 
_ ae 
ratio —" for 
€o 
smaller 


The nonspherical particles is 


calculated but 
here the relation is obeyed quite well. 


somewhat than even 
In the 
derivation of vaN HEERDEN’s equation it is 
shown that the shape factor of irregular particles 
approaches 1 if effective particle diameter is 
used instead of average particle diameter and it 
has been tacitly assumed above that the same 
RvuBey’'s 
of the experimental and 


holds for equation. A comparison 
‘aleulated vaiues of 
Re, for nonspherical particles, shows this assump- 
tion to be incorrect, whereas it is found to be 
reasonably accurate in the case of Re,. Therefore, 
the experimental ratio probably deviates from 
the calculated one because of the absence of a 


shape factor in Rusey’s equation. Despite these 
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Fic. 4. 


XY 


deviations the fit of the experimental points 
to the calculated line is sufficiently good to give 
confidence in the validity of equation (5) for 
fluidization numbers above = 3 x 10*. Whereas 
Rusey's equation is obeyed for lower fluidiza- 
tion numbers, VAN HEERDEN’s equation is not, 
as has been pointed out above, and the experi- 
mental points suggest that the true curve should 
the indicated by dashes which 


follow line 


Re 
T for very small 
€9 

very large particles. 


approaches a constant ratio 


fluidization numbers, i.e. 


If equation (5) were correct in that region a 


Re 
point would be reached where —' would become 


€¢ 

less than 1, ie. particles would be elutriated 
without fluidization having set in. The values 
of Witnetm and Kwavk also lie close to the 
modified curve and demonstrate the validity of 
the curve for fluidization both in gas and in liquid. 

It is interesting to the physical 
significance of the curve (Fig. 5). Its plot shows 
five distinct branches. The first branch for very 
large fluidization numbers virtually 
constant ratio, due to the dependence of both 
Re, and Re; on the cube of particle diameter. 
This is the region of Stroxes’s law and the ratio 


consider 


indicates 


Microspherical alumina 
@ From data of Witnetm and Kwauk 
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Microspherical catalyst 










































































Rey 
0 
law (equation 1) is divided by equation 4. This 
branch will be limited by the validity of Srokes’s 
law. (Srokes’s law is not obeyed by particles 
of less than 1 micron diameter). The ratio of 76 
represents the greatest possible flexibility between 
critical and terminal Reynolds numbers (and 
therefore mass velocities) attainable in a fluidized 
reactor without overhead loss of particles. 
Section 2 of the curve is the region of change of 
terminal velocity from turbulent 
with a gradual change of slope, while Re, remains 
proportional to 4. 


= 76 can also be obtained if Sroxkes’s 


viscous to 
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Critical and terminal velocities in fluidization 


Eventually the terminal velocity becomes fully 
turbulent and Re, will then be proportional 
to d*. If the limitations on the van HEERDEN 
equation previously referred to did not exist, 
this would result in the straight line portion 
of slope 1-5 starting around section 3. In actual 
fact section 3 of the curve is very short and may 
be only an inflection point; the slope then 
changes again in section 4 as the critical fluidiza- 
tion velocity changes from viscous to turbulent 





10 
80 

















10 
Porticle diameter 
Fic. 6. Effect of particle density. 
Fluidizing gas-nitrogen. Temperature 20°C. 


with a change of the proportionality of Re, 
from d;3 to d/*. When the critical velocity 
is fully turbulent both Re, and Re, are propor- 
tional to d,'* and their ratio consequently has a 
constant value of about 10, as in part 5 of the 
curve. 

It is of interest to note that van HeEeRpEeN 
et al [5] proposed values of 70 for the streamlined 
region (corresponding to Section 1 of the curve) 
and 8-5 for the fully turbulent region (i.e. Section 
5); these are in good agreement with the present 
results. 


The relationship of Ber may have practical 
Reo 

implications in the design of a fluidized catalyst 

reaction vessel, where elutriation of particles is to 

be avoided and it may also be useful in the des gn 


of filtration equipment in fluid catalyst reactors 
with particles of a wide range of diameters. 
Its practical implications may be discussed 
briefly. Flexibility, that is the number of times 
the critical fluidization velocity may be exceeded, 
varies from 76 for the smallest to about 10 
for the largest particles. The actual diameters, 
for which these values are attained, depend on 


solid density and the density and viscosity of the 
fluidizing medium. To reach the value of 76, 
smaller particles are required for heavier catalysts. 
(See Fig. 6). Fluid density has a similar and 
fluid viscosity an opposite effect, hence increasing 
temperature with gas or vapour as the fluidizing 
medium causes the fluidization number to increase 
rapidly for any given particle diameter, i.e. 
greater operating flexibility is possible at higher 

10 

80 















































Fic. 7. Effect of fluidizing gas temperature. 
Particles — Glass ‘ Ballotini.’ Gas-Nitrogen. 


temperatures (see Figure 7). Flexibility is also 
increased with spherical catalyst compared with 
ground catalyst of the same effective diameter. 
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NOTATION 
diameter of particle m 
effective diameter m 
arithmetical mean diameter (of narrow diameter 
range) m 
2 
ae 
Po (Ps — po) 4° 8 
acceleration due to gravity 9-8 m/sec? 
Go = critical fluidization velocity kg/m? sec 
Gy =terminal mass velocity kg/m? sec 
1 
QF 
Go d G 
=~! = 7 of for spherical particles 
M M 
: Gr d, Gr d 


fluidization number 


for spherical particles 
“ “ 
=linear terminal velocity m/sec 
dynamic viscosity kg/m sec 
density of fluidizing medium kg/m* 
= density of solid kg/m* 
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Circulation in liquid drops 


(A heat-transfer study) 


P. H. Catperspank* and I. J. O. Korcurnskit 


Chemical Engineering Department, University of Toronto 
( Received 4 June 1956) 


Abstract—The process of material and heat transfer to and from liquid drops moving through 
another continuous fluid has been examined. Continuous phase heat-transfer coefficients have 
been measured for the system mercury-aqueous glycerol solutions in which the mercury con- 
stituted the disperse phase. These data are in good agreement with those previously reported 
for solid spheres up to a Reynolds number of about 200, at which point drop oscillation commences. 


The resistance to heat transfer of the disperse phase, expressed as an effective thermal diffusivity, 


VOL. has also been measured, for the system bromobenzene-aqueous glycerol solutions. These results 
6 were found to be in general agreement with much reported work and with the mathematical 
956/57 model of Kronic and Brink which applies to liquid drops undergoing internal circulation. 


Some observations concerning the rigid, circulating, and oscillating behaviour of drops are 
made. The conclusion is reached on both experimental and theoretical grounds that circulation 
within a liquid drop results in an effective mass or thermal diffusivity of 2-25 times the molecular 
value. This result applies for the whole régime of circulation which exists between rigid sphere 
and oscillating behaviour. The onset of drop oscillation results in a considerable increase in 
this factor. Some effective diffusivities are reported for the oscillating régime which has not 
however been sufficiently studied to enable firm conclusions to be drawn. 


Drag coefficients for falling liquid drops in the above-mentioned systems are also reported 
and reveal the limitations of recent correlations in this field. 


Résumé—Les auteurs examinent le mécanisme du transfert de matiére et de chaleur vers ou a 
partir de gouttes liquides se déplagant dans un autre fluide continu. Les coefficients de transfert 
thermique dans la phase continue ont été mesurés pour le systéme mercure-solution aqueuse de 
glycérine dans lequel le mercure constitue la phase dispersée. Les valeurs obtenues concordent 
avec celles déja obtenues pour des sphéres solides jusqu’éa un nombre de Reynoups d’environ 
200, valeur & partir de laquelle commence loscillation de la goutte. 

Exprimée sous forme d'une diffusivité thermique effective, la résistance & l'échange thermique 
de la phase dispersée a été également mesurée pour le systéme bromobenzéne-solutions aqueuses 
de glycérine. Les résultats sont en concordance satisfaisante avec la majorité des travaux décrits 
et avec le modéle mathématique de Kronic et Brink qui s'applique 4 des gouttes liquides 
présentant une circulation interne. 

Les auteurs ont fait quelques observations relatives au comportement de gouttes, rigides, en 
circulation ou en oscillation, A partir de résultats expérimentaux et théoriques, ils arrivent a la 
conclusion que la circulation a l'intérieur d'une goutte liquide entraine une diffusivité thermique 
ou massique d’environ 2,25 fois la diffusivité moléculaire. Ce résultat s‘applique a la totalité du 
régime de circulation s‘étendant de la sphére rigide a la sphére oscillante. L’amorgage d’oscillation 
de la goutte entraine un accroissement considérable de ce facteur. Les auteurs donnent quelques 
valeurs de diffusivité effective pour un régime oscillant qui n'a pas encore été suffisamment étudié 
pour permettre des conclusions. définitives. 

Les coefficients de trainée pour des gouttes liquides tombant dans les deux syst¢mes mentionnés 
ci-dessus indiquent les limitations des corrélations récentes dans ce domaine. 


* Present address : Chemical Research Laboratory, DSIR, Teddington, Middlesex, England. 
Tt Present address : Shell Development Company, Emeryville, California, USA. 
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INTRODUCTION 


A riGoRous analysis leading to the prediction of 
the performance of spray columns in any of their 
varied applications requires a knowledge of the 
following quantities : 


Drop size. 

Drop velocity. 

Continuous phase heat or mass transfer 
coefficient. 


Effective thermal or mass diffusivity 
within the drop. 


Drop size data have been correlated by Hey- 
worth and Treypa. [1] in terms of nozzle 
diameter and the physical and flow characteristics 
of the system. Drop velocities have been corre- 
lated most recently by Hu and Kinrner [2]. 

Continuous phase heat and mass transfer 
coefficients for systems involving solid spheres 
have been reported and are summarized by 
McApbams [3] and Suerwoop and Picrorp [4}. 

If the drop is assumed to behave as a rigid 
sphere in which transfer rates are solely deter- 
mined by the value of the molecular diffusivity, 
it becomes possible to design a spray column from 
the foregoing the Grober 
chart [5]. This chart presents in convenient form 


the solution to the mathematical problem of 


correlations using 


the continuous and disperse phase 


A restriction imposed by the use of 


summing 

resistances. 
the Grober chart is that the temperature or con- 
centration obtaining in the continuous phase be 
constant. For most practical cases, this condition 
would not be realized. If, however, the equili- 
brium and operating lines are linear, the solution 
of AMUNDSEN [6] may be used. The latter method 


is prohibitively laborious and a convenient general 


solution perhaps employing stepwise methods 
would be most desirable. However, even if such 
a method were developed, a serious objection 
remaining would be caused by the observation of 
Garner [7] and others, that the liquid drop is 
generally in a state of internal motion. The use 
of the molecular diffusivity instead of some 
enhanced effective value, would 
expected to lead to serious overdesign. 

The purpose of the present paper was that of 


therefore be 


evaluating the continuous phase-resistance and 
the effective diffusivity obtaining during the 
unhindered fall of drops. 

In the circulating régime, it was demonstrated 
that the continuous-phase transfer resistance 
follows closely that reported for solid spheres, 
while the effective diffusivity within the drop is 
a constant multiple of the molecular diffusivity. 


General Principles 
} 


It was decided, at the outset of this work, to study 
the process of heat rather than mass transfer 
within drops as previous workers have clearly 
demonstrated the difficulties in studying the pro- 
cess of mass transfer where trace contaminants, 
by their accumulation at the interface, lead to 
inconsistent results. Every precaution was how- 
ever taken to ensure scrupulous cleanliness of the 
fluids and equipment. 

Drop sizes and terminal velocities were meas- 
ured directly as was the continuous phase-heat- 
transfer coetlicient. The continuous phase-heat- 
transfer coeflicient was evaluated in a series of 
experiments in which mercury drops were allowed 
to fall through hotter or colder aqueous solutions 
of glycerol. The resistance to heat transfer of the 
mercury drops was calculated to be a negligible 
fraction of the total resistance in these experiments 
in which a wide range of Reynolds and Prandtl 
numbers was obtained by variation of the amount 
of glycerol in the solution. It was found that the 
heat transfer during drop formation could be 
practically eliminated by allowing the drops to 
form in air at one drop diameter above the con- 
tinuous phase. 

In a further series of otherwise identical experi- 
ments, bromobenzene was used as the disperse 
phase. In this case the resistance to heat transfer 
of the drops was large. The resistance to heat 
transfer of the disperse phase was deduced from 
the measured overall resistance by making proper 
allowance for that of the continuous phase, as 
evaluated in the experiments with mercury drops. 
As it did not prove feasible to form the drops 
above the continuous phase, a correction had to 
be applied to allow for heat transfer during drop 


formation. 
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Circulation in liquid drops 


Drag Coefficients 
Drag coeflicients were determined for mercury 
and bromobenzene drops falling through aqueous 
glycerol solutions. 

Drop velocities were then deduced in the usual 
manner from a plot of f- (Re)? against Re where 
4 gd (pa — Pe) Pe. 


3 yu? 
In the experiments with mercury drops, a 


for spherical drops /f-: (Re)* 


number of nozzles were calibrated so that the 
drop size obtained in the heat-transfer experiments 
was known. When bromobenzene was used as the 
disperse phase, it was found more reliable to 
measure the drop size during the course of the 
heat-transfer experiments. 


The Continuous- Phase Resistance 


It may be easily shown from a heat balance, 
integrated for the case of a constant continuous 
phase-temperature, that the continuous-phase 
heat-transfer coefficient may be deduced from 
experimental data by means of the relationship : 


Cpa pa’ d vg — %, Ti — To 


h . 
c 6 "+. a6 


(1) 

Further, h, may be correlated in the usual way 
with the physical and flow properties of the 
system from the functional relationship : 


Nu = f (Re, Pr) 


Disperse-Phase Resistance 


Under conditions of constant continuous-phase 
temperature and negligible continuous-phase 
resistance, the equation for diffusion in a rigid 
sphere, due to NewMan [8] is : 


Ti— To _ 
Ti —Te 


VERMEULEN [9] showed that this result could be 
accurately represented by the empirical equation: 


antt\4 
gBalt—e or) (4) 


Kronic and Brink [10] included the effects of 
internal circulation on diffusion and for Re < 1 
and at zero-interfacial tension derived : 


gpa1-* Pee (5) 
8 nel 
where the eigenvalues A, and yp, have been 
tabulated [11]. 

DanckweEnrts [12] has pointed out that equation 
(5) may not be as strictly limited as the above 
restrictions suggest. 

Equation (5), which applies under conditions 
of constant continuous phase temperature and 
negligible continuous phase resistance was found 
in the present work to be well represented by the 
empirical equation : 


2-25 2 ny 


E= (1 ~~ = (6) 


This fact is demonstrated in Fig. 1 in which 
equations (3) and (4) and equations (5) and (6) 
are compared. 
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Equations for diffusion into stagnant and circu- 
lating drops. 
2 
t 
E vs. *", PLOT 
r 
d 
equation (5); © equation (4); 
\, equation (6). 


equation (3); — 


If equations (4) and (6) are compared it may 
be readily seen that internal circulation of the 
type envisaged by Kronic and Brink may be 
represented by a constant effective diffusivity 
equal to 2-25 times the molecular value. 


Finite Resistance in Both Phases 


Groeer [5] has solved the equation for diffusion 
into a sphere in the case where the continuous 
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phase-resistance is finite. These results have been 

put into a convenient form and plotted in Fig. 2. 06 
The restriction of constant continuous phase- 04 

conditions still applies to the proper use of Fig. 2. 
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Fic. 2. Equations for diffusion into spheres with finite 
continuous-phase resistance. O02 


Transfer during Drop Formation 


The extent of heat transfer during drop formation 
may be determined approximately by plotting 








experimental values of In(l — E,) against 





the column length for experiments in which the ) 
latter is varied, and extrapolating the curves to 
zero column length. ‘) Plot of — log (1 — Ep) vs. L for bromobenzene. 
Drops falling in 74:3°% glycerol solution, 
10 Ep = 0-12 
] ] | @ Nozzle No. SS-1. Nozzle No. SS-2. 
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O04 


3. Heat transfer during bromobenzene drop formation 
as determined by extrapolation procedure. 


02 — t ’ This procedure is not rigorous since a linear 
driving potential is assumed but was adopted in 
this work in default of a better method. 
The fractional approach to equilibrium during 
drop fall may be found from : 
| Ey, —_ E, 
“M si J 
1— Ep 
0 t 


cm 

















. — As previously stated, it was found possible to 
Plot of — log (1 i) vs. L for bromobenzene. f . : 
~~ S ae form mercury drops just above the continuous 
Drops falling in 56-5% glycerol solution. ’ F 
Ep = 015 phase-surface and to thereby reduce the experi- 
@ Nozzle No. SS-1. © Nozzle No. 22-8. mental value of E, to zero. 
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In the case of bromobenzene, the above pro- 
cedure led to distortion as the drop moved 
through the continuous phase-surface and drops 
had therefore to be formed below the surface when 
finite values of E, lying between 12 and 16% 
were recorded (see Fig. 3). 


EXPERIMENTAL DETAILS 


Physical Properties of Liquids 
Reagent ACS grade 97-8%, glycerol by weight was 
made up into solutions of various concentrations 
as determined by refractive index measurements. 
The specific heats of the solutions were obtained 
from Perry [13], the densities from the inter- 
national critical tables [14], the thermal con- 
ductivities from the data of Bares [15] and the 
VOL. viscosities from the data of Suee.ry [16). 
6 Mercury was taken from the supplier and 
1956/57 redistilled. Physical properties were obtained 
from the Handbook of Chemistry and Physics {17}. 
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(a) Interfacial tension between mercury and glycerol- 
water solutions. 
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(b) Interfacial tension between bromobenzene and glycerol 
water solutions. 


Fic. 4. Interfacial tension between mercury and glycerol 
solutions and between bromobenzene and glycerol solutions. 


Interfacial tensions between mercury and glycerol 
water solutions were determined using the drop- 
weight method of Harkins and Humpurey [18] 
results being shown in Fig. 4. 

Bromobenzene was obtained from BDH (boiling 
range 154°C to 158°C). The refractive index as 
measured was 1-5599 compared with the value of 
155977 quoted in the International Critical 
Tables [14]. 

Interfacial tensions of bromobenzene in aqueous 
glycerol solutions were measured by the pendant 
drop-method [19]. The solubility of bromobenzene 
in glycerol solutions was found to be too small 
to be detected. 


Nozzles 

Five glass nozzles were prepared from capillary 
tubing and precision ground to a flat face, the tip 
nozzle diameter being measured with a toolmaker’s 
microscope. Two sharp-edged stainless-steel 
nozzles were also used. The internal diameter of 
the nozzles ranged from 0-0328 cm to 0-4320 em. 


Drop Sizes 

Drop sizes were measured by observing the volume 
of disperse phase formed when a measured number 
of drops had coalesced. As drop rates varied from 
1 to 5 per second, a photocell and electronic 
counter were used to count the number of drops. 
Mercury drops were obtained varying in size from 
about 0-2 to 03cm diameter. Bromobenzene 
drops lay in the size range 0-386 to 0-540 cm 
diameter. 


Drop Velocities 

Drop velocities were determined by photographic 
means. By the use of a “stroboloom” in a 
darkened room, several photographic exposures 
at regular time intervals of a drop falling through 
the continuous phase were obtained. 

The frame of the photograph also included an 
accurate linear scale. Results, plotted in terms of 
drag coefficients are shown in Fig. 5. 

High-speed photographs were taken with a 
“ Fastax’ camera and showed that terminal 
velocities were attained at a point very close to 
the nozzle. It was also seen that mercury drops, 
formed one drop diameter above the continuous 
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(a) Friction factor plot for mercury drops falling through 
glycerol-water solutions. 
© exp. data. ——— solid spheres line. 
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(b) Friction factor plot for bromobenzene drops falling 
through glycerol-water solutions. 
© exp. data. ——— solid spheres line. 


Fic. 5. Drag coefficients for mercury and bromobenzenc 
drops compared with the behaviour of solid spheres. 


phase were not distorted in passing through the 
liquid surface. Bromobenzene drops were however 
appreciably distorted and this fact led to the 
decision to form them below the continuous 


phase-surface. 


Fic 6. 


Apparatus for heat-transfer 


measurements. 
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(a) Heat transfer apparatus. 


Elevation section 
——_——— 6cm ————_——o~ 


— 42cm ——~ 


| .— 3°Sem _,., 
0:42cm 





A, 


ht. Drop 
catcher 



































Stopper 








Polythene tubing 
Exit disperse phase 
thermistor 
Continues phase 


To burette 
To reservoir 


Column wall 





Disperse 
phase 
outlet 


Inlet 
ees 

\ phase 
Plan . thermistor 


(b) Column bottom assembly. 
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Fic. 6 (continued) 
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(c) Column top assembly. 


The effect of wall proximity on the terminal 


drop velocity was also studied by using columns 


of varying diameter. It was found that the effect 
of column diamcter had disappeared when the 
column diameter was about ten times that of the 
drop. 

A column diameter of 3-71 cm was used in the 
heat-transfer work and the effect of wall proximity 
could be discounted. 


Heat-transfer Measurement 


A schematic diagram of the equipment used is 
shown in Fig. 6. 

A vacuum-jacketed glass column was fed with 
the disperse phase from a constant head stainless- 
steel tank. The disperse phase was collected in an 
hourglass-shaped catcher at the base of the 
column and coalesced in the narrowest section 
where a thermistor was located. The level in the 


Constant 
head 
bubbler 





-Vent 








+ Coil 


+ Thermometer 


Filling line 





.e— Tank wall 
































—Nipple 








— 
(d) Disperse phase feed tank. 


catcher was maintained constant by means of a 
“ swinging-leg ” interface control. The disperse- 
phase temperature in the feed tank was measured 
with a calibrated thermometer reading to 0-02°C. 

The continuous phase reservoir consisted of a 
thermostat bath controlliy z to + 0-02°C over the 
range 15°C to 30°C. The continuous phase was 
pumped from the thermostat in a loop through a 
calibrated rotameter, upwards through the column 
and back to the thermostat. 

Temperatures were read by means of thermistors 
at the points indicated in Fig. 6. The thermistors 
could be switched in turn into a resistance net- 
work and enabled temperatures to be read with a 
+ 0-02°C. The thermistors and ther- 
mometers were calibrated against a National 
Research Council of Canada standard thermo- 
meter. The thermistor bridge circuit is shown in 


Fig. 7. Thermistors were obtained from the 


precision of 
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Fic. 7. Bridge circuit for use with thermistor temperature- 
measuring elements. 


Victory Engineering Corporation, New Jersey 
(Veco type 51-Al resistance at 25°C ~ 10° ohms) 
and showed excellent stability after pretreatment 
at 100°C for one month. 

The effective column length was 20-5cm in 
most of the runs. Temperature readings were 
taken after about half an hour of steady operation, 
the rotameter reading noted and drop sizes and 
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(a) Hu-Kintner correlation for mercury drops in glycerol 
water solutions : 
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Hu-KIntNer data. © Present work. 


velocities checked. The difference between the 
temperatures of the in- and out-flowing continuous 
phase was in no case greater than 0-03°C so that 
the continuous phase could be assumed to be at a 
constant temperature without sensible error. Care 
was taken to use a column length which would 
ensure an accurately measurable-temperature 
difference between the phases at the base of the 
column. 


RESULTS AND Discussion 
Drag Coefficients 
Experimental drag coefficients have been reported 
in Fig. 5 and show a deviation from solid-sphere 
behaviour only at Reynolds numbers above 200, 
when drop oscillation was observed to occur. 

A visual study of the behaviour of bentonite 
clay particles suspended in the bromobenzene 
drops clearly revealed the fact that internal 
circulation was fully developed at Reynolds 
numbers well below 200, so that the effect of 
internal circulation on the drag coefficient, was 
too small to be detected. 

A comparison between the present results and 
the correlation suggested by Hu and Kintner [2] 
is made in Fig. 8. The latter workers established 
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(b) Hvu-Kinrner correlation for bromobenzene drops in 
glycerol-water. 
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Fic. 8. Comparison between experimental results and the correlation of Hu and KintNer [2]. 
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Circulation in liquid drops 


a functional relationship between the drag co- 
efficient and the Reynolds and Weber groups by 
dimensional analysis, and determined the ex- 
ponents of the above groups by a trial and error 
fit of experimental data. It was thus found 
possible to correlate the data for falling liquid 
drops in water by means of a single line plotted 
on the co-ordinates shown in Fig. 8. From a 
comparison with the present work it is clear that 
the effect of continuous phase viscosity, not 
investigated by Hu and KiInrNer, is considerable 
and can lead to discrepancies when the field liquid 
is not pure water. 


Continuous- Phase Heat-Transfer Coefficients 


Values of the continuous-phase heat-transfer 
coefficient were calculated from the data obtained 
with mercury drops by means of equation (1). The 
continuous-phase temperature was taken as the 
average, the maximum change being 0-03°C. The 
heat loss from the mercury drop to the air during 
its formation was calculated to be insignificant. 
The heat loss at the coalescing section was not 
known but assumed to be insignificant as the 
contact area and driving force were both small. 
The continuous-phase velocity was taken as the 
average velocity of the fluid in the column. A\l- 
though the continuous phase flow was laminar, 
a very flat velocity profile appeared to exist when 
the notion of suspended particles was observed 
under strong illumination. This appears to result 
fortuitously from the design of the entry and exit 
sections. | Continuous-phase heat-transfer co- 
efficients are plotted in Fig. 9 and compared with 
published results for heat transfer to solid spheres 
and the analogous situation of mass transfer to 
liquid drops where the dispersed-phase resistance 
is negligible. The range of the variables covered 
in the present work was : 

Re — 05 to 3,000 

Pr — 6 to 3,300 

Nu — 10 to 120 


It may be seen that the transfer coeflicients for 
liquid drops depart from those found for solid 
spheres at a Reynolds number of about 200. As 
this point coincided with observed drop oscillation 
it is apparent that oscillation causes some bound- 


ary-layer degradation and an increased transfer 
coeflicient. 
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Fic. 9. Continuous phase heat transfer coefficient correla- 
tion from mercury drop data. . 
Exp. data. 
Curve from McApams [3]. 
Data of TANG, et al. [28]. 
X Data of Licur and Pansine [21]. 
@ Data of Coutson and SKINNER [20]. 


Disperse-Phase Heat-Transfer Resistance 


Experiments were performed using bromobenzene 
as the disperse phase under conditions which 
ensured that a large part of the overall resistance 
to heat transfer lay in the bromobenzene phase. 
Corrections for heat transfer during drop 
formation were made as already explained (Fig. 3). 
This correction was found to be equivalent to a 
fractional approach to equilibrium due to drop 
formation of between 0-12 and 0-16 and is of 
the same order of magnitude as the equivalent 
results of Licut and Conway [22] and HEEerTJEs 
et al. [11] which were between 0-18 and 0-20. 
Values of the continuous-phase heat-transfer 
coefficient (h,) were deduced from the results 
obtained using mercury drops (Fig. 9). From the 


parameter > and the experimental E,, 


Pa: Pa: 

values, the effective thermal conductivity of the 
drops was found by means of the modified 
Groper chart shown in Fig. 2. The ratio of 
effective to molecular diffusivity (R) was then 
calculated. 

A range of Reynolds numbers of from 10 to 
600 was covered, results being shown in Table 1, 
together with those calculated from data in the 


literature. 





Table 1. 
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1. J. O. Korcuinski 


R values calculated from this work and from literature data. 





Source of data 


Garwin and Smrra [23] 
Garwin and SMirn [23) 
McDowe .t and Myers [24) 
McDowe.t and Myers [24] 


This work 


Heertyes et al. {11} 


This work 


Hucues and GiLuicanp [25] 


GaRNER and SKELLAND [26] 


- 


PANsinG [27 


| Cold benzene drops rising in 


System 


hot-water heat transfer 

Hot benzene drops rising in 
cold-water heat transfer 

Cold kerosene drops rising in 
hot-water heat transfer 

Cold xylene drops rising in hot- 
water heat transfer 
Bromobenzene drops falling in 
hot or cold glycerol-water- 
solutions heat transfer 
Isobutanol rising in 
water-mass transfer of water 


drops 


R range 


Re range Type of drop behaviour 


| 


1-72-2-1 1500-1780 No drop oscillation 
730-810 
circa 
2000* 
circa 
2400* 


1-6-2-4 No drop oscillation 
No drop oscillation 


No drop oscillation 


10-120 | No drop oscillation 


20-200 No drop oscillation 


Bromobenzene drops falling in | 


hot or cold glycerol-water 


solutions heat transfer 


315-620 | Drop oscillations present 


Water drops falling through | 


carbon dioxide atmosphere- 


mass transfer of carbon di- | 


oxide 

Nitrobenzene drops contain- 
ing acetic acid falling in water- 
mass transfer of acetic acid 
Perchloroethylene drops con- 
taining acetic acid falling in 
water-mass transfer of acetic 
acid 


_ Drop oscillations present 
Type of drop behaviour not 
stated by authors. Oscillations 
indicated by R values 
Drop oscillations present, al- 
though oscillation decay indi- 
cated by range of R for con- 
stant conditions 


300-3000 
270 
370 
480 
375 
595 


700 





* Exact continuous-phase temperature not given by authors ; 


reference [24] was used to estimate the Reynolds number. 


In the circulating régime R values found in 
this study varied from 1-8 to 3-3 a result in general 
agreement with corresponding figures calculated 
from literature data and supporting the value of 
R = 2-25 deduced from the Kronic and Brink 
model. 

In the oscillating régime, R values found here 
and deduced from published data are considerably 
higher than is found for circulation alone. 

The constancy of R with column length, which 
was found, appears to disprove the hypothesis of 
oscillation decay advanced by Hucues and 
GILLILAND [25] and gives strong support for the 
fundamental soundness of the use of an effective 
diffusivity to express the transfer resistance of 
the disperse phase. 

A serious limitation of the present work is its 
failure to define the range of conditions over 


temperature of 170°F given in sample calculation of 


which circulation is present and the point at which 
oscillation starts. Much more experimental work 
is needed before this issue can be clarified. 


CONCLUSIONS 


The influence of internal circulation in liquid 
drops is to raise the effective diffusivity to a 
value equal to 2-25 times the molecular diffusivity. 

No precise data is yet available to define the 
régime over which pure circulation occurs. 

Oscillation causes the effective diffusivity to be 
greatly increased above the value obtained for 
pure circulation. 

Existing correlations for the drag coefficients 
of liquid drops require to be extended to account 
for wider variations in the viscosity of the field 
liquid, 

Internal circulation by itself has a small, and 
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in this work unmeasurable, effect on the drag latter effect is doubtless the manner in which 

coefficient or the continuous-phase heat-transfer ‘ pulsed’ columns achieve a superior perform- 

coefficient, the latter two quantities agreeing with ance. 

solid-sphere behaviour. Acknowledgements—One of us (1.J.0.K.) wishes 
Drop oscillation increases the drag coefficient, to thank the National Research Council of Canada 

the continuous phase heat-transfer coefficient and for financial assistance in the form of two student- 

the effective diffusivity within the drop. The _ ships. 


NOTATION 
Ayn = an eigenvalue. 


= specific heat of drop, cal/g°C. 
drop diameter, cm. 
= mass diffusivity, em? /sec. 
= base natural logarithms, dimensionless. 
fractional approach to equilibrium, dimensionless. 


fractional approach to equilibrium, during drop formation, dimensionless. 


VOL, ; fractional approach to equilibrium during drop fall, dimensionless. 


total fractional approach to equilibrium, dimensionless. 
9 56/57 = drag coefficient, dimensionless. 
gravitational acceleration, 980 cm /sec?. 
continuous phase heat transfer coefficient, cal /em? sec’C. 
= thermal conductivity, calem /cm*sec°C. 
column length, em. 
= an integer, 1, 2, 3... 00, dimensionless. 
Nusselt number, hed /k, dimensionless. 
= Prandtl] number, Cpy/k, dimensionless. 
drop radius, cm. 
diffusivity or conductivity ratio, dimensionless. 
Reynolds number, dvp/yu, dimensionless. 
Schmidt number, »/ Dp, dimensionless. 
Sherwood number, kd/D, dimensionless. 
contact time, sec. 
continuous phase temperature, °C. 
average temperature, of drop, °C. 
average temperature of incoming drop, °C. 
average temperature of outgoing drop, °C. 
continuous phase velocity, cm/sec. 
relative velocity between drop and continuous phase cm/sec. 
= thermal diffusivity, em? /sec. 
effective thermal diffusivity, em*/sec. 
= 3-142. 
= continuous phase density, g/cem*. 
pq = disperse phase density, g/em*. 
te = continuous phase viscosity, poise. 


#, = an eigenvalue. 
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APPENDIX 


Continuous-Phase heat-transfer coefficient correlation. 
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Calculations for the bromobenzene-drops — glycerol-water-solutions heat-transfer experiments. 





h, 
Run He __ 
No. | poise f (Re? Re em? °C sec 
| 





| 0-0258 | 8-02 x 10*| 338 0-0993 . 0-365 0-524 0-433 
0-0258 8-21 338 0-104 5- | 0476 | 0-614 0-559 0-525 
0-0260 | 8-65 350 0-105 | 15- 1-04 1-33 0-762 0-283 
0-0270 | 7-19 315 0-105 i 2-94 3-89 0-969 0-037 
0-0272 | 8-30 340~— |S 0-0960 i | 2-89 3-30 0-956 0-052 
0-0272 8-33 340 | 0-102 3 4-28 5-17 0-976 0-029 
0-0222 | 10-4 385 0-113 3B | 416 5-97 0-977 0-028 
0-0222 | 19-0 | 530 «| O116 | 0-313 | 0-377 0-476 | 0-624 
0-0192 | 26-1 620 0-121 9: 0-605 | 0-752 0-641 | 0-427 
0-0260 | 15-0 465 0-110 | . 0-891 0-986 0-720 0-333 
0-0270 | 12-8 : | 0-105 | 19 2: 2-69 | 0-908 0-109 
0-0270 | 14-7 d 0-103 “ | 2 | 2+ 0-893 0-127 
0-0274 | 144 4 | 0-103 9! | 3 3-7: 0-953 0-056 
00-0274 13-5 | 0-102 | “4 | 3-06 ° 0-938 0-073 
| 0-102 0-0503 3 296 0-390 0-718 
0-102 0-0509 2- 983 | 0-55: 0-496 0-593 
0-0991 0-0487 2- BE ‘7 0-551 528 
0-100 0-0500 | 128 3-9: 2- | OT749 295 
0-100 0-0498 ‘ 2: 0-831 199 
‘101 00-0508 2- Or 3-3 0-869 “154 
101 +6 . 0-0472 2: Bf | 3 0-849 178 
100 f 118 | 60-0526 : : WM 0-336 781 
100 | 106 0-0537—s| | 0-8: zo 0-429 672 
100 . 107 0-0541 3: 22 | O6 0-521 563 
1000 | «124 112 0-0532 3° 3. 69 | O716 334 
100 ‘27 » 115 | ©0529 5 3 | 4-5 0-662 398 
100 “2! 114 | 00530 +60 | { 2-3: 0-776 263 
100 3 118 0-0527 5S | 2-3 0-753 291 
(0215 | 1-93 > 31-5 0-0336 | 9%: 0-653 | 0213 | 0-327 ‘765 
“414 f 11-3 0-0272 1-54 . 0-467 606 
424 3s 11-8 0-0300 2 0-573 485 29 
435 0-543 12-0 0-0274 5-5 1:3 | 0-838 184 87 
1 
1 


























a & 
2H Ae 
see 


> =) by 


te 


59 
‘79 


“45 


™ 


te tw te 


te 


435 0-545 | 120 0-0271 0-830 | 0-193 83 
“436 0-425 10-2 0-0283 0-910 0-102 | ? 1-94 
0-439 0-484 10-8 0-0271 O-916 0-096 r+ 2°31 
“437 0-807 16-3 00-0264 | | 0-165 0-312 0-782 4-6 2-34 
“437 0-788 16-1 0-0273 } OK “lk 0-250 0-367 0-720 4+: 2-37 
438 0-805 16-3 | 06-0268 88 | 21: 0-455 | 0-465 | 0-608 5 2-06 
0-434 0-858 I73 | 0-0262 09 | 5-96 1-25 0-693 | 0-349 5 | 241 
| 0-433 0-915 | 181 | 0-0265 *! 5-6 114 | 0668 | O377 | 15 2-33 
0-435 | O-782 | 160 | 00264 | 6 . 2-07 | 0-789 0-240 . 2-08 
0-476 | 0-917 18-2 0-0266 | 1-77 | 0-831 0-192 . 3-05 





Nore: Ey values for the various continuous phases were : 28-7% G Ey 0-160; 565% G- E, = 0150; 743 % 
G — Ey = 0-120. 
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Distillation at minimum variable reflux 
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Abstract— Minimum heat input in a distillation in which the reflux is variable is calculated 
by a modification of UNprerwoop’s minimum-reflux calculation. The method is based upon the 
analytical theory of distillation previously given by one of the authors. Instead of the molal 
liquid and vapour rates actually occurring in the column, generalized rates are defined which are 
dependent on the molal heat rates of the vapour phase and of the liquid phase. The equations 
between the generalized rates are formally identical with the equations between the molal rates 
in a distillation at constant reflux. To obtain these generalized rates the molal vapour heat 
content must be approximated by a linear function of the equilibrium constant of a reference 
component K and the molal liquid heat content must be approximated by a linear function of 


K These approximations need, however, only be valid in the relatively narrow temperature 

VOL. , 1 aliifisae attus walk —— re 
range between the two “ pinches.” Constant relative volatilities have also to be assumed for 

6 that region only. 

Résumé—Les auteurs présentent une méthode pour calculer la quantité minimum de chaleur 

nécessaire pour effectuer une distillation a reflux variable. La méthode est fondée tantét sur 

la méthode d°'UNpDERWwoop pour calculer le reflux minimum en cas d'un reflux constant, tantot 


sur la théorie analytique de la distillation développée antérieurement par un des auteurs. Au 


956/57 


lieu des courants de vapeur et de liquide actuels on introduit des courtants généralisés qui dépend- 
ent des enthalpies de la vapeur et du liquide dans la colonne. Ces quantités généralisées sont 
liées les unes aux atures par des équations identiques formellement aux équations gouvernant 
la distillation a reflux constant. Pour obtenir ces quantités généralisées il faut admettres que 
lenthalpie de la vapeur puisse étre représentée par une foncti n linéaire de la constante d’équilibre 
K dun des composants du mélange, pris comme substance de référence ; il faut également que 


lenthalpie du liquide puisse étre représentée par une fonction linéaire de I] suflit que ces 


approximations soient valables entre les deux régions de composition invariable (pinches). II 


faut aussi que les volatilités des composants aient un rapport constant dans cette intervalle de 


température seulement 


OssecT or THE PaPer case of a large difference in temperature between 
UnprRrwoonp’s well-known method for the caleula- top and bottom of a column, 


tion of the minimum reflux ratio for distillations 


UxnpERWoopD’s METHOD 


in which the components have a constant relative 


volatility and in which the reflux remains constant 
throughout the column is widely used at present 
[1, 3). 
that the method can be generalized to calculate 


It does not seem to be generally known 
minimum reflux in a distillation in which the 
components still have a constant relative volatility 
but the reflux is variable. 

In this paper the theoretical basis of the 
generalized method is given, and the theory is 
illustrated by an example. The discussion also 
explains why the assumption of a constant relative 
volatility still gives satisfactory results even in 
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We may recall that in the simple case in which the 
distillation is carried out in such a way that only 
two components are present in specified quantities 
in both products, the distillate and the bottom 
product, the method consists in the calculation 
of a root common to two algebraic equations 
one pertaining to the top section, the other to 
the stripping section. The equations are* 


+ HY; 
Z; Yip _ R, n (1) 


a— w 


*A list of notations is given at the end of this paper. 
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for the top section and 


2p 
———"--* 
for the stripping section. 

Instead of using concentrations the 
unknown w, the present authors prefer the use 
of molal rates and of the variable ¢ = w™ as 
many of the calculations become simpler in that 
case. Equations 1 and 2 rewritten in molal 
rates and introducing ¢ as the unknown variable 
become after some rearrangement : 


and 


(4) 


The condition of minimum reflux is that 


(3) and (4) have one common root which is thus 


eqs. 
also a root of the equation : 
(5) 


This root is the one situated between the values 
x,' and a,' if h and | are used to denote the two 
distillate and in the bottom product. The index 
h refers to the heavy key, the index / to the 


components " that are specified in the 


light key a, %,). 
The difference L, L 
content of the feed; it is equal to the liquid 


depends on the heat 


portion of the feed under feed tray conditions*. 

Let t 
the 
reflux corresponds to the values of L, or of L, 


2 be the root of eq. (5) which satisfies 


condition ay" < 2.< a", then minimum 


as calculated from eq. (3) and eq. (4) respectively 
with f Q. 


Variable reflux. 

It has the 
authors that the calculation of any distillation 
in which the reflux is variable but the relative 


been shown by one of present 


volatilities are constant, can be reduced to a 


calculation of a distillation with constant reflux 


*The difference L, — L, is often calculated in a different 
ic. from a heat balance over the column assuming 


way, 
constant reflux in both sections. 


and constant «,;, provided that the molal heat 
content of the vapour and of the liquid can be 
approximated by a certain quadratic function 
of the equilibrium constant of the reference 
component [4, 5]. The reference component is 
the one for which « = 1. Van Wisk and pr Haas 
vAN Dorsser have discussed some specializa- 
tions of the general expressions [6]. The proof 
that a distillation variable reflux 
calculated as a distillation with constant reflux 


with can be 
will now be repeated here for a special approxi- 
mation of the 
given in the last mentioned article. This approxi- 
mation is used in the following example of 


molal heat content which was 


calculation of minimum variable reflux. 
Let us assume that the heat content per mole of 
liquid H, and per mole of vapour H, in the 


column are given by the expressions 
H, Hy HY K (6) 
i, HS Hq. K (7) 


in which A is the equilibrium constant of the 
reference component and the H’s in the right- 
hand sides are constants. 

In a distillation with variable reflux the 
vapour and liquid streams have to satisfy the 
heat balance, which for the section B indicated 


1 has the form 


H, = H,.V.—H 


om m im + 


in Fig. 


_L (8) 


’ 


Here V,, is the total molal vapour flow leaving 


the mth tray so V,, 2.0 f Uw 


the molal vapour flow of component i. Similarly 


Ln 1 
Substituting eqs. (6) and (7) in eq. (8) we 
obtain among others a term A, V,, and a term 


L.,,/K 


m m+tl 


denotes 


, 
po ae 


The former is equal to 2, V, ,,/x; 
the latter is equal to 2,L,; ,,,,%,- In fact 
denoting the absorption factor of the reference 
component on the mth tray by A,, = L,,/V,, K,, 
and its reciprocal value by S,, (stripping factor) 
and remembering that the equilibrium relation 
on an ideal tray demands that L, ,, = 27" A,, Vj. 


we have 
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nN=WN 
n=N+*1 m=Mt1 














Fic. 1. Diagram showing the numbering of the trays in 
the two sections of a distillation column. 


A y; aS 


m+i~i 


- 
2a, L 


m+ L, m+ 
i.m+l* 

Thus eq. (8) can be written as : 

H, = HSV, + AE” ZV; n/%: 

AOL, — A LD, wi % (9) 
in which K,, or K,,,, do not occur. Either 
one of the molal flows V, ,, and L,; ,,, , can be 
expressed in the other one. 

Eliminating for example V,,, by application 
of the relation L; ,,,, = V;,, + B,; one obtains 
Hy, = 2;([H,% — Hf + A/a, — Af] - 

_ z, [H HS x, B,. (10) 


A, = [H, — HY 


HS /a; — H{” «,] 
constant (11) 
we can rewrite (10) as 
Ly = 2A, L,_4., = A, 4 
+ &; [AO + A x,| B; 


constant. (12) 


From the mathematical point of view eq. (12) 
is the condition for distillation at constant molal 
reflux of a mixture with “ generalized ” liquid 


flow rates, which are defined as 


L en (13) 


Introducing the generalized residue rates 
By = A, B; (14) 
and the generalized vapour flow rates 


\ (15) 


(o,m) 


The equations expressing the conversation of 
matter 


B. (16) 
hold also for the generalized components or 


I » = V, B,) (17) 


“(i,m + (i,m) 


as follows from eq. (16) by multiplication of both 
sides with the constant A;. Eq. (12) is equivalent 
to eq. (8) which expresses the conservation of heat. 
Similarly from the equilibrium conditions 


L i ie (18) 


‘1m 


we obtain the generalized equilibrium conditions 
on multiplication by A, 


L be, Vera (19) 


(i,m) 


From eqs. (12), (17) and (19) it follows that a 
distillation at variable reflux can be treated 
mathematically as a _ distillation at constant 
reflux in which the generalized flow rates occur 
instead of the actual flow rates. The equation 
of conservation of matter and the equilibrium 
conditions, in which the actual relative volatilities, 
the actual absorption factors, and the actual 
stripping factors occur, are exactly equal to the 
corresponding relations in the actual distillation 
but the equation for the heat balance is based 
upon the validity of the expressions (6) and (7) 
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Table 1. Separation Scheme of Sizx-component system*. 





liq. feed 


0-250 
0-078 
0-180 
0-095 
0-038 
0-025 


0-010 
0-012 
0-070 
0-075 
0-072 
0-095 


Methane 
Ethane 
Propane 
Butane 
Pentane 
Hexane 
0-666 


Total 0-334 


vap. feed 


F; D, 


0-260 
0-090 
0-246 
0-003 


20-6 
5-07 
2-06 
1-00 
0-429 
0-206 


0-260 
0-090 
0-250 
0-170 
0-110 
0-120 


0-004 
0-167 
0-110 
0-120 
0-401 


1-000 0-509 





* The relative volatilities are average values estimated from different sources for a temperature of 83°C and a 


pressure of 20 atm. 


for the molal heat content. If these expressions 
are valid the actual absorption 
calculated from the distillation 
generalized reflux is constantt. 


factors are 


in which the 
A similar conclusion holds for the top section. 
The heat balance applied to the section A in 
Fig. 1 gives 


Hy as +1 Vi (20) 


| 
We can derive from eq. (20) 
J pare Hy 


lin 
5, (Hi + HY? a) D, 


r 


() 1 


(21) 


constant 


(generalized equation of conservation of heat). 
In addition we have the equations 


Ve D L 


lisa +d (a) lin) 


(generalized equations of conservation of matter) 


Vin) x, 5, Lin) (23) 


(generalized equilibrium conditions). 


Consequently eqs. (3), (4) and (5), which result 
from eqs. (16) and (18) together with the con- 
dition of reflux 2, L, ,, = constant, 
are automatically also valid in terms of the 
15, 17 12). 


constant 


generalized quantities (eqs. and 


tit is obviously incorrect to calculate the equilibrium 
constant of the reference component from the absorption 
factors inserting the generalized flow rates since by 
definition A,, = ZL; ,,/K 2a; 3 and this is in general 
not equal to 21); 4) /K Lat for the same value 
of K. 


(i,m) 
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of 


has 


generalized 
a previous 


The theory distillation 
quantities 
communication to which reference is made for 


further details [4, 5]. 


using 


been given in 


CALCULATION OF MINIUMUM 


CONSTANT REFLUX 


EXAMPLE OF 


The theory will now be applied to a six com- 
ponent system which was previously discussed 
by Jenny [2]. Pressure 20 atmospheres, tempera- 
ture of the feed 81°C. The component rates are 
in moles for one mole of feed. 
light key component, butane is the heavy key ; 


Propane is the 


butane is chosen as the reference component. 
L, + 0-334 
and eq. (5) possesses the root 2 = 0-719334 
which is in the range a," = 0-485 to ag! = 1-00. 
This value of 2 inserted in eq. (3) gives the 
0-553 and L, 
reflux conditions. These liquid rates correspond 
0-553 0-599 0-923 
the bottom 


For constant reflux we have L, 


results L, 0-887 for minimum 
to a top reflux ratio R, 

vapour reflux 
1-212. 


and ratio at 


Ry 


a 


Tue 
The molal heat content of the liquid and of 
the vapour in the column must now be known. 


MoLaAL ENTHALPIES 


For the calculation of the present example 
a constant specific heat of 0-6 keal kg has been 
assumed for the liquid. Thus if the liquid 
enthalpy is arbitrarily put equal to zero at 
t = 0°C and M is the molal weight, 
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H, = 0-6tM. (24) 


The molal enthalpy for the vapour is assumed to 
be 


i, (0-5t + 70) M. 


Both expressions are in kcal kg mole and ¢ is 
in degrees centigrade*). 

The relation between the molecular weight 
of the two phases and the temperature or the 
equilibrium constant of the reference component 
must be obtained from a correlation. As such a 
correlation was not available for the present 
case the following procedure was used. The 
molecular weight is known at the top and at 
the bottom of the column; it is also known 
how the M of pure aliphatic hydrocarbons 
varies with their boiling temperatures under 20 
atmospheres. The latter values were plotted as 

*As the object of this paper is to explain the method 
for calculation of a distillation with variable reflux rather 
than to calculate a given separation rigorously, the 
authors have not attempted to calculate the most accurate 
values for the heat contents involved. 





Pure 
cliphatic — 
hydroc 


Molecular weight —.+ 








05 TO 15 “3-0 


arm 





Fic. 2. Molecular weight as a function of the equilibrium 

constant of the reference component in the distillation 

of the six-component system. These molecular weights 

are used to convert enthalpy per unit of mass to molal 
enthalpy. 


a function of the temperature and this graph 
was used to construct two similar curves almost 
parallel to the first one, on which the points 
corresponding to the terminal conditions are 
situated. The results of this procedure are 
given in Table 2 and are shown in graphical 
form in Fig. 3 as a function of K and of 1/K 


S000 








§ 
8 





kcal/kg mole —+ 
8 
8 

















Vk — 
Fic, 3. Liquid molal heat content as a function of the 
reciprocal value of the equilibrium constant of the reference 
component. The straight line is the approximation of H, 
between the two pinches used in the text. The dotted 
vertical lines indicate the values at the pinches and of the 
feed. 


Table 2. Molal enthalpies of the liquid and of the 
vapour in the column as a function of the tempera- 


ture or of K. (kcal/kg mole). 





K m.w.liq. I, 
0-255 41 885-6 30 
0-30 1084 
0-40 1529 
0-50 1958 
0-60 2363 
0-70 2767 
0-80 3168 
0-90 3557 

1-00 3965 
1-10 4355 
1-20 4763 
le 5187 
l- 
l- 


m.w.vap. 
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5586 


6594 
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respectively. The molal heat of the vapour can 
excellently be represented by the expression of 
eq. (7), but the molal heat of the liquid must be 
approximated by different expressions of the type 
of eq. (6), each valid in a restricted range of 
temperatures only. 

REFLUX 


MINIMUM VARIABLE 


Unperwoop’s relation is based upon the fact 
that an infinite number of trays is present in 
each section and that it is still possible to have 


one common tray in both sections which is the 
feed tray [3]. 
streams V; , in the top and L, , in the bottom 


The composition of the terminal 


is of no importance for the condition of minimum 
reflux and neither is the composition of the reflux. 
Since an infinite number of trays is present 
between the regions of invariant composition 
(pinches) and the feed tray in each section, it is 
sufficient to use eq. (6) for the molal liquid heat 
in the the 
pinches. Although the absorption factor of the 


region of K values between two 


reference component is approximately known 
beforehand in both pinches, the K values can 
only be calculated from them after the value 
found.* 


one has to select a certain interval of A values 


of minimum reflux has been Thus 
more or less arbitrarily, and adapt the constants 
in eq. (6). Afterwards it must be verified whether 
the K values at the pinches are actually situated 
in that interval. 

A value K 0-63 or 1 A 


to the degree of vaporization of the feed. The 


1-59 resp ynds 


straight line in Fig. 3 has been drawn such that 


*In the theory set forth in ref. (5), it has been shown 
that in the stripping section of a column the product 
of absorption factors up to and including that of the 
mth tray, Ap (m), is equal to the coefficient of ? in the 
development of the following expression into a power 
series in ¢ 
2; Vicig/A x,t) 


A 
Lo) — 2; 8/0 — =f) 


' (t) 
The product of absorption factors is defined as 


Ap (m) y 2 Ay: ---- - 


The above expression can be expanded into elementary 
rational fractions, so that one may write 


A 


in which ¢,~' is a root of eq. (4) (for generalized quantities) 

and ¢, is the residue of the function A if ¢ approaches 
¢; ! or 

lim A (1 

t=@;-! 


¢; ¢). 


Several conventional methods exist for the calculation of 
the residues, but we shall not discuss them now. The 


absorption factor product at the mth tray is thus 


Ap (m) = 


2; C; ¢;”. 


The absorption factor A,, = Ap(m)/Ap(m — 1) accord- 
ing to the definition of Ap (m). 

Since in the column section below the pinches (inclusive), 
the unspecified components do not appear, eq. (4) does not 
contain the components 1 and 2 and the maximum «a 
occurring in that equation is 2, = 2-06. The largest ¢ 
found from that restricted equation is obviously equal to 
the absorption factor of the reference component in the 
lower pinch, since it is the limiting value of 

2; C; ¢;™ 
‘ 2; Cj $;™ . 
This ¢ depends on the value of Ry. However, in the 
which, as 


A ~?¢ 


max fOr m approaching infinity. 


present case, By, is a small fraction of 2; By, 
can be shown from a discussion of eq. (4), has the conse- 
quence that one root is close to a, '. This is the largest 
root and therefore the absorption factor in the pinch is 
close to ag, since ¢ is equal to the reciprocal value of 
the smallest root. 


In the top section the theory shows that the function : 


Ss 1 


max 


Apay'tt Ap (2y' @ ... Ap (ny? it" 


takes the place of the function A in the stripping section, 
and one obtains the analytical expression 
Ss a 


(t) 


“11, 0) ql 
2; Diy (1 


Or expanding into elementary rational fractions 


a 
‘a — 4,0 


from which follows that 
Ap (ny Zc; (6;')". 


The ¢,s are the reciprocal values of the roots of the 
denominator in the first analytical expression of S and 
(¢;y'. 

If the unspecified components in the top section 5 and 
6 are ignored the largest value of ¢;' corresponds to the 
smallest root of the restricted equation. This root is 
somewhat smaller than 1-00 in the present case, owing 
to the fact that D,,, is a small fraction of 2; D,;,. Since 
one has in the pinch of the top section 


the c,"'s are the residues of that expression for ¢ 


A 1 
lim - pi 
n= 


Ss _ 
ao Ap(n) 


pinch ~ ? max 


the stripping factor in this pinch is somewhat higher than 


one. 
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3. Constants for calculating generalized rates. 





A; 


35670 9274 
6509 — 585-8 -—— 
2139 8-556 526-2 15-42 


2640 140-9 
8116 892:8 
20152 2418 


3760-3 


A; is defined by eq. (11), B,;, and Dj) 


| 
(1, + H/o] By) [H,{© + H, «,) D; 





— 8833 
— 366.9 
+ 422-4 

7-920 1068 + 11-27 

1427 

3053 


10394 5563 


by eq. (14) and as A, D; respectively. 





it differs less than 10 per cent from the molal 
liquid heat in the range of 1 A values 0-95 to 
2-05. The constants in the right-hand sides of 
eqs. (6) and (7) have been calculated from this 
line for the liquid and from the corresponding 
graph for the vapour, giving the equations 


H, = 5683 — 1925 K (26) 


H 1458 1940 K 


t 
The resulting coefficients A, are given in Table 3 


together with other data necessary for the 
calculation of minimum reflux. 
Returning now to eq. (5) we insert 


Fi = By 


L,—L the expression 


Li) Li Li) Vin i (i) 
Hy + 2; [H° HY 2; | dD, 
+ 2, [H+ HL /a,) B; + 2, Dy 


D,,, instead of F,; and instead of 


ZL, Din = He 


The constants H, and Hp are the net upwards 
heat flows in the stripping section and in the top 
The difference H, — Hy 
is, ignoring heat loss, equal to minus the heat 
0-334 « 2480 
Thus 


section respectively. 
introduced per mole of feed, or 
— 0-666 x 4660 3932 k cal kg mole. 


Ly) —_ Ly = 
= $250. 


~ 89382 — + + 103894 


The root of eq. (5) which lies between a,' = 1 
and a;' = 0-485 has a value = 2 = 0-684912. 
If this root is inserted in eq. (4) in which the 
generalized residue rates B,, are substituted 


5458 is found at minimum reflux 
(y — 5458 — 8258 = 2200. 


The net upward heat fluxes in the column are 


a value L,,) 
and from this follows L 


now found from eqs. (12) and (21). One obtains 
2 D,,: 


recalling that V,,, = L,,) 


H, = 5458 — 5563 105 (29) 
H ,, = 2200 + 10394 — 8767 + 8827 (30) 
These values represent the net upward heat 
They 


characterize the separation just so as the value 


flux at minimum variable _ reflux. 
of minimum reflux characterizes a_ separation 
in which the reflux is constant. They are 
constant in the entire column and also valid 
outside the range of validity of eq. (26). Thus 
H,, is the difference between the heat content 
in the vapour returned to the column from the 
reboiler V, and the heat content of the liquid 
flowing down from the lowest tray, Z,. Similarly 
H,, is’ the difference between the heat content 
of the vapour from the upper tray and that of 
the liquid reflux at the top of the column. 
Rertux Ratio’s at Cotumn TERMINALS 

The reflux ratio’s at the top and at the bottom 
are now calculated assuming equilibrium conden- 
sation and vaporization respectively. The com- 
positions of the liquid reflux from the condenser 
and of the vapour reflux from the reboiler are 
given in Table 4. 

For the determination of the reflux ratio R, 
the molal vapour reflux V, must be known. 
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Table 4. Reflux compositions for equilibrium in the 
reboiler and in the condenser. 





7 


Dj/x; YiB #ip 
0-0126 
0-0177 
0-1190 
0-0030 


0-082 
0-116 
0-782 
0-020 


0-0082 
0-1670 
0-0472 
0-0247 


0-033 
0-676 
0-191 
0-100 


0-2471 


0-1523 1-000 1-000 


The concentration of comp. i in the vapour reflux is 
Yip = % B K/g49, with K the equilibrium constant of 
the reference component in the reboiler; as DT y;, = 1-00 
one has K = 1-62. Similarly z;,) = D;/a; K * (-599 and 
K = 0-255 in the condenser. 





It is calculated by solving eqs. (16) and (8) 
substituting HW, = — 105 subject to the condition 
that K, at the lowest tray (m 
such a value that the liquid from this tray is 
at its boiling point or 


1) must have 


2; [vy %: Vin B, x, | Ky Vy 
The set of values K, = 1-34 and V, = 0-517 is 
found to satisfy these conditions, resulting in 


Ry, = 0-517 0-401 1-29 (31) 


This is nearly equal to R, for constant reflux. 

In the top section, three similar conditions 
3827, 
eq. (22) for ordinary i.e. nongeneralized flowrates 


must be satisfied viz. eq. (20) with H/,, 


and the condition that V, is at dew temperature. 
One obtains A (for n 1) = 0-0875 
Lo = 0-750 resulting in 


and 


R,, = 0-750 0-599 = 1-26. 


(32) 
This is a considerable difference with R, = 0-923 
for constant reflux. 

In a distillation with variable reflux the reflux 
ratio itself is, not a 
quantity ; instead of it enter the net heat fluxes 
in both sections in the column, which are coupled 
through the heat input in the feed and through 
possibly occurring heat losses. 


however, characteristic 


Eeuitiprium Constant at THE “ Pincn” 


The temperatures in the regions of invariant 
composition, are now calculated. 
It is necessary that both temperatures fall within 
the range of applicability of eq. (26). In the 
region of invariant composition the temperature 
remains constant, and thus the vapour heat 
and the liquid heat in eq. (8) or in eq. (20) have 
to be taken at equal values of K in each equation. 
The stripping factor of the reference component 
in the pinch of the stripping section is approxi- 


* pinches,” 


mately equal to the smallest root in eq. (4) 
components and 
The stripping factor in 


ignoring unspecified using 
generalized flow rates. 
the pinch of the top section is equal to the 
greatest root occurring in eq. (3), if again un- 
specified components are ignored and generalized 
Thus one has exactly 


two equations to determine the total liquid or 


quantities are inserted. 


vapour rate and the equilibrium constant of the 
For 
instance in the top section one obtains with 
L,,, = 2200 a value 1-0072 for the greatest root 
of eq. (3), thus KV /L = 1-0072 and eq. (20) 
becomes 


reference component for each “ pinch.” 


8827 = [H,(K) — H,(K)|V 4 


H,(K) »« 0-599. (33) 


The two equations are solved by trial and error. 
A value of K is assumed and the molal heats 
are calculated from the expressions (26) and 
(27). The vapour rate V resulting from eq. (33) 
is then inserted in the stripping factor together 
with L V — 0-599. One obtains K = 0-536 
or LK 1-867 and V = 1-28. By the corres- 
ponding procedure for the stripping section one 
obtains K 0-820, or I1/K 1-22 
V = 0-567. The stripping factor of the reference 


and 


component in the pinch of the stripping section 
is the smallest root in eq. (4) inserting L,,, = 5458 ; 
it is 0-489. 

The values of K corresponding to the invariant 
compositions are apparently both well within 
the range of validity of eq. (26), (cf. Fig. 3). 

Since the calculation of minimum heat flux 
is based upon the fact that and infinite number 
of trays between the two pinches is present, 
the relative volatilities of the components in 
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the range of temperatures between the two 
pinches are only essential. A large variation of 


a; 


outside this region does not 


affect the 


conclusions at all. 
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Ap (m) 


i 
H,, ne 


H 
HL, HS 1) 


NOTATION 


function of a variable ¢ defined by the 
relation 

A 1 + Ap(ijt 
Ap(m)t™ +.... 
absorption factor of the reference 


Ap (2) 


component A on tray number 
m. 

absorption factor of component i on 
tray number m. 


absorption factor product for the 


L 
AV 


reference component at tray number m 
Ap (m) = A, A,. ove Me 

relative volatility of component i 
in respect to the reference component 
a, = K,/K. 

residue rate of component i. 
generalized residue rate of component 
i, defined as B,» = A; B;. 

subscript referring to the bottom 
section. 

constant occurring in the expression of 
the absorption factor product 

Ap (m) Z, ¢; $;". 

distillate rate of component i. 


_generalized distillate rate of component 


i, defined as D,,, = A; Dj. 

feed of component i. F; = D; + B,. 
reciprocal root of denominator in the 
analytic expression of A. One has 
Ap (m) 2; ¢; ¢;"- 

net upward heat flux in bottom section. 
net upward heat flux in top section. 
molal liquid heat content. 

constants occurring in the expression 
of molal liquid heat content as a linear 


1 
functio T= 
unction 0 K 


molal vapour heat content. 

constants occurring in the expression 
of molal vapour heat content as a 
linear function of K. 

subscript to indicate heavy key. 


L 


J J 


(i,m) " (1, ”) 


Vin 


um’ 


(i, m)° Luin) : 


subscript to indicate a component. 


= equilibrium constant of the reference 


component. 


= value of K on tray m. 


equilibrium constant of component i 
on tray m. 


= total liquid rate in bottom section. 


total generalized liquid rate in bottom 
section. 


- liquid rate of component i from tray 


m or n respectively. 
generalized liquid rate of component i 
from tray m or n respectively. 


= total liquid rate from tray m or n 


respectively. 

total liquid rate in top section. 

total generalized liquid rate in top 
section. 

subscript indicating light key. 

factor defined by the relation. 
HL HH, HS» a; — HH; a;. 
The component rates are transformed 
into generalized rates by multiplication 
by A 


molecular weight. 


= index for tray number in bottom 


section. 

index for tray number in top section: 
reflux ratio in condenser. 

vapour reflux ratio in reboiler. 
stripping factor of the reference 
component on tray m or n respectively. 


= sommation over all components. 


mathematical variable used in the 
calculation of the absorption factor 
product. 

temperature. 

total vapour rate from tray m or n 
respectively. 

vapour rate of component ¢ from tray 
m or n respectively. 

generalized vapour rate of component 
i from tray m or n respectively. 
generalized total vapour rate in top 
section. 

matheinatical variable occurring in the 
equations for calculation of minimum 
reflux. 

molal concentration of component i 
in liquid residue. 

molal concentration of component i in 
gaseous distillate. 

characteristic root occurring in the 
equations for calculation of minimum 
reflux. 
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Boundary conditions of flow reactor 
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Abstract—An analysis is presented of the boundary conditions for a steady-state flow reactor 
with axial diffusion and first order reaction. Conclusions regarding reactor properties are reached 
as the result of simultaneous solution of the three differential equations for the reaction section, 
fore and after sections. Axial diffusion may occur in the latter two sections. The limits of 
behaviour, represented by a plug-flow reactor without diffusivity, and a fully-stirred reactor, 
are dissussed. Contributions of the diffusivities in the three parts of the system to the course 
of reaction in space are considered. 


Résumé— Les auteurs présentent une analyse des conditions limites pour un réacteur a écoule- 

VOL. ment stationnaire avec diffusion axiale et une réaction de premier ordre. Leurs conclusions 
6 concernant les propriétés de réacteur sont le résultat de la solution simultanée de trois équations 
956/57 différentielles concernant la section principale de réaction, la section d’entrée et la section de 


sortie. La diffusion axiale peut intervenir aussi bien dans ces deux sections que dans la section 
principale. Ils discutent les limites de comportement correspondant d'une part aux réacteurs 
& piston sans diffusivité et d’autre part aux réacteurs a agitation parfaite. Ils considérent 
l'influence propre sur le progrés de la réaction dans l’espace, de la diffusion dans les trois sections 


du réacteur. 


Tue analysis of a plug flow chemical reactor 
which includes disappearance of a substance 
by chemical reaction and its transport in the 
axial direction by means of diffusional processes 
is of interest. The diffusional mechanisms may 
be molecular or turbulent. In the limit of in- 
finitely large diffusivity the solution of this 
problem is identical with that of a fully-agitated 
chemical reactor. With zero diffusivity, the 
solution approaches that for a plug flow reactor, 
the type of reactor frequently assumed to be 
valid in purely chemical studies. 

Two previous papers dealing with the problem 
are by Danckwerts [1] and HuLBurt [2]. Present 
work discusses mathematical conclusions that 
result from the choice of different sets of boundary 
conditions. 

The steady state differential equation in 
dimensionless form is : 


1 df(z)_ df (2) py 
ca 2° " 


L , 
where the Peclet group, Pe > L is a length 


term that characterizes the system, here taken 
for convenience as the bed length. Differently 
defined Peclet groups may be introduced if 
desired, using appropriate ratios of characteristic 
length terms. AR is the rate group for a first 


. kL , ; 
order reaction —, and / is the fraction of reactant 
u 


remaining. It is assumed that there is no volume 
change during reaction. z is a dimensionless 
distance in the axial direction varying between 
z=0 and z= 1. 

DaNnckKwenrts [1] used the boundary conditions: 


1 df 
J(0—) Pe dz 


b. (1) = 0 (3) 


where {(0 —) in equation (2) represents the 
fraction remaining in the feed stream as it 
approaches the bed entrance, /(0 +), just 


f(0 +) (0+) (2) 


*Present address: Chemical Engineering Department, The Catholic University of America, Washington, D.C. 
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within the bed entrance. The difference between 
these two values arises because of diffusion within 
the bed, represented by the third term. There 
is a discontinuity at this boundary as reported 
which will be shown later to be not necessarily 
correct. The second boundary condition was 
established on intuitive grounds. 

Hu.surt (2) employed the following set of 
boundary conditions : 


a. f(0—-)=f(0+)=1 (4) 
b. Same as equation (3). 


Equation (4) neglects the diffusion directly 
within the reactor. 

The solutions resulting from the choice of 
boundary conditions between these two authors 
differ by a constant factor f (0) defined in equation 
(22) of this paper. 

In considering a flow reactor it seems desirable 
not to neglect the role of the fore and after 
sections. In the present model these two end 
sections (Fig. 1) stretch from © to + oO. 
The reactor bed is situated between 0 to 1 on 
the This method 


dimensionless = co-ordinate. 











of treatment introduces boundary conditions 


at 2z + oc which seem to be intuitively 


acceptable. In turn, the boundary conditions 
at z = 0, 1 are clarified by means of the behaviour 
of the solutions of the differential equations in 
the fore and after sections. 

The differential equations for the three sections 


are : 
1 af df 
Pe,d2 dz 
1 df df 


Pe, d* dz 


1 af df 
Pe.d2 dz 


0 


0 


The Pe groups may differ between sections 


, eye ae 
f= N, + Nz exp Pez 


Pe 
{=i = 


for reasons of velocity, presence of particles 
and the like. 
The six required boundary conditions are : 


(a) t( 
u, A, |r(o )— Pe, Zio ] 
1 df 


Pe, dz 


20) 1 


uy, Ay | f(0 +) 


The total continuity equation states : 


u, A 


a a 


u, A, = u, A,. (10) 


The use of the continuity condition permits 
simplification of the boundary condition for 
section a, equation (9) becoming : 


1 df(0 —) 
b) 0 : : 
(D) f' Pe, dz 


ja+- 22 


(13) 


(e) f(l —) (14) 


(f) f() is finite. (15) 


Conditions (b) and (d) result from the con- 
servation of reactants at exit and entrance of 
the bed, taking into account flow and diffusion. 

Conditions (c) and (e) arise from the intuitive 
argument that at steady state the concentration 
should be continuous between sections at the bed 
entrance and exit. 

The general solutions for the three sections 
containing the arbitrary constants, N;, and 
obtained by solving equations (5), (6) and (7) 
are : 


(16) 


‘ Pe 
ta)z+N exp 4 
2 


N, exp Pez 
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After inserting boundary conditions the com- 
plete solutions for the previous equations are : 


1-f 

1 —f (0) 

ee, Perr) | (a Pe, (1 — 2) 

f go exp (“5 i ' a) exp ( —a =} 


a Pe,(z —1))\} 
=) 


= exp (Pe,z) z <0 (19) 


(1 — a)exp( 0<z<1 (20) 

and 

f= f (1) = 2ag, exp (Pe,/2 

where, a = Vl + 4R/ Pe, 
2 


=constant,z >1 (21) 


£o 


—,— 
aPe, 


—— = 
(1 + a)texp( 2%) — (1 - ajPexp (— = 


and 
J (0) 
: Io Ip 
£o ‘a { a)exp(“—*) (1- a)exp| — - (22) 


In comparing the solutions for the three sections 
of the bed it is interesting to observe first that 
the f — z profile in the fore section (a) depends 
upon the (a)—(b) boundary concentration, 
f (0), which in turn depends in its absolute value 
upon the parameters in the reaction section, 
Pe, and R. However, the profile in the centre 
reaction section (see equation 20) does not depend 
upon conditions in the fore section (a). The 
flat profile in the third section (equation 21) 
is the result of setting constant N, in the general 
equation (18) equal to zero in order to satisfy 
the boundary eondition in equation (15). The 
level of concentration, f(z > 1) in this section 
depends upon reaction and diffusion in the centre 
section. The value of the diffusivity in the after 
section does not enter into the results. 

It now becomes of interest to discuss the 
boundary between sections (a) and (b). The 
following are pertinent points : 

1. The value of the flux in the fore section 
is constant at a value of 1-0 which may be shown 
by substituting into the expression for the total 
flux, f(z) — c. df (2) 
Pe, dz 
equation (19). This finding remains true at the 
boundary, equation (11). 


.the solution for the section. 


Increasing 
Pe, 
Section Section 


(a) (b) 





2. For f(z< 0) the curve has the general 
shape as shown in Fig. 2. Thus as the Peclet 
group is caused to assume very large values, 
approaching conditions of zero diffusivity, the 
profile in section (a) approaches a step function 
(Fig. 2). Pertinent compositions and ‘gradients at 
f(0) are indicated in Fig. 3. 


Section 


The solution for the reaction section 
(equation 20) is identical with that of 
Danckwerts. Both Danxckwerts and_ the 
present authors used a value of unity for the 
flux at the inlet boundary. However in the 
present instance this value was _ established 
from an analysis of the fore section. The boundary 
condition as written by Danckwerts is the 
equivalent of the step function mentioned above 
1 df(z< 0) 
Pe, 


The boundary conditions between the reaction 


since he dij not include the term 


and after section also were the same in the present 
work as that of Danckwerts. 

Let us consider now the argument by which a 
zero gradient at the end of the reaction section 
is the logical limit. The solution for section ¢, 
cquation (21) yielded a flat profile. Thus, in the 
toundary condition equation (13) the last term 
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is zero. Using also equation (14), it follows 
GF(l'-)_»o 
dz 


immediately that In effect, the 


boundary condition at z= + o is drawn to 
the end of the reaction bed. 

An additional feature of interest regarding the 
gradient at z = 1 is the rationalization between 


df (1 —) 


~ 


the above result, namely that 0, 


and the conclusion that stems from consideration 
of a reaction with plug flow but without diffusion. 
In this latter case the differential equation is as 
follows : 


df — Rf 


| (23) 


From the well-known solution of this equation, 


Equation (24) for reaction without diffusion 
arises as Pe, + © in equation (25). This limit 
must be taken within the bed, 0< z< 1. Then 
asz Il, df > 0. 

dz 

For any finite value of diffusivity, ie., Pe, 

df 
" dz 
development is {1 — exp (a Pe, (z — 1)}. 

Fig. 4 portrays concentration profiles through 
the three sections of the reactor for selected 
values of Pe and R. The logical limits of plug 
flow reaction without diffusion (Pe, = o) and 
of infinite diffusivity as with a fully ‘stirred 
reactor (P, = 0) are shown in section b. We 
note especially the finite gradient for the Pe, = oo 


< ©, asz>1 0. The key term in this 








Pe,= 2-667 


+ + 


Pe,= co 


Section c 


Section b 

















+0 -O08 -O6 -0O4 


2 


Zz 


Fic. 4 


f may achieve a value of zero only at infinite bed 
length. At the finite bed length, z = 1, 
l 

o R exp (— R) 


This apparent conflict in boundary gradients for 
and without may be 
resolved by examination of the derivative of 


(24) 


solutions with diffusion 


equation (20), which is as follows : 


df 
dz 


$R(exp 


(1 +a)? —(a- 


Pe, (1 —a)z 
af a) )fa exp(aPe, (z 


1)? exp (— a Pe,) 


curve at z= 1. For any finite diffusivity the 
gradient at that position becomes zero. The fore 
section (a) was arbitrarily assumed to have 
Pe, = 1 in the illustration and the profiles as 
shown result. Because the Peclet groups were 
not taken with equal values between the sections 
the slopes are discontinuous at the boundary. 
The Hvucsurrt solution may be arrived at by 
normalization through division of _ the 
Danckwerts solution by /(0). Although con- 
venient for estimating relative importance of 
reaction and diffusion rates the solution is not 
consistent with the limit of a stirred reactor 
with infinite diffusivity in the reaction section. 
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APPENDIX* 


Boundary conditions (11) and (13) are, in general, 
crue for a three-dimensional and non-steady-state reactor 
since they contain terms which represent the only 
nechanisms considered for flow of reactant across a 
youndary. 

When no diffusion occurs in the fore section (Pe, > «) 
voundary condition (11) becomes : 


ld 
Yo 4) (26) 


J -) =f +) - 5. gt 


This equation holds during a time-varying change in f 

and the resulting discontinuity does not disappear in the 

steady state as evidenced by the present solution. 
For boundary condition (13) : 


fQ—-)=f(+) 
if 


and, as Pe, > @ 


in the nonsteady state if Pe, is finite. This result may be 
shown as follows : 


yO)... » Wm <e 


As (Pe),.> @, = = 


This means physically, however, that section c acts on 
the surface z = 1 as a “ sink,” a result which is clearly 
absurd since no reaction takes place in section c. It 
Y=) _ 

_ 2 


follows therefore that as (Pe). ~, 0 


of(l +) 


although, in general — + 0. Thus, even in the 


1 
dz 


nonsteady state, no discontinuity will develop for any 
period of time at boundary (1). At the entrance boundary 
(0), however, the discontinuity is set up and maintained 
by the reaction and diffusion within the bed. If no 
diffusion occurs within the bed, no discontinuity will 
occur at the entrance boundary, (0). 

It is to be noted that equation (27) is valid only for 
Pe, = © in the non-steady state, but for all values of 
Pe, in the steady state. Also, Pe, must be finite for 
equation (27) to be valid in either case. 


*The generality of these equations and the two special cases were pointed out to the authors by Prof. ANDREAS 
Acrivos of the Department of Chemistry and Chemical Engineering, University of California, Berkeley 4, 


California. 


NOTATION 


@ = a group, (1 + 4R/Pe,)'? 
A = cross-sectional area through which flow occurs. 
D = diffusion constant in Fickian equation, molecular 
or turbulent. 
Jf = fraction of reactant remaining 
J (0 —) = composition at upstream side of fore section- 
reactor boundary 
{(0 +) = composition at downstream side of fore 
section — reactor boundary 
f(1 —) = composition at upstream side of reactor — 
after section boundary 
Jf (1 +) = composition at downstream side of reactor - 
after section boundary 


is first-order reaction rate constant 
L = reaction section bed length 
N = constant in general solution for differential 
equation 
Pe = Peclet group, Lu/D 
R = reaction group, kL/u 
u = linear velocity 
z = dimensionless distance in the axial direction. 
Normalizing constant is L. 
Subscripts 
a refers to fore section 
b refers to reaction section 
ec refers to after section 
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Book reviews 


Principles of Chemical Engineering Thermo- 
dynamics. Ernest D. Witson and Haroun C. Ries : 
McGraw-Hill, 1956. pp. ix. + 376. 56/6s. 


Tue authors describe their book as being written for an 
undergraduate course in thermodynamics for chemical 
engineers. The contents include a discussion of the first 
and second laws, with applications to fluid flow, steam 
power plants, refrigeration, phase equilibrium and reaction 
equilibrium. Graphical representations of thermodynamic 
data are included, together with numerous useful illustra- 
tions in the text and about 120 problems for the student 
to solve for himself. 


The first six chapters, which are concerned with pro- 
perties of matter and the applications of the first law, 
are taken quite slowly and at a rather elementary level. 
The second law is thus not reached until almost half-way 
through the text. This results in a much faster pace in 
the second half and the last two chapters in particular 
are too highly condensed: phase equilibria, ideal and 
non-ideal solutions, reaction equilibrium and the third 
law, are all compressed into the space of 59 pages. These 
important topics are thus given a rather summary treat- 
ment. 


The authors make very extensive use of the perfect 
gas law in the development of the theory and in their 
preface they remark that they believe this to be the best 
method of achieving simplicity of presentation. (And 
of course many other textbooks follow the same pattern.) 
Whilst admiring the simplicity which has been achieved, 
the reviewer feels that it is a mistaken policy in 
thermodynamics teaching to give undue prominence to 
perfect gases ; its effect is that students, having become 
habituated to the use of certain equations applicable only 
to perfect gases, blithely apply them to any type of sub- 
stance, liquid, solid or gaseous. Thus students commonly 
write dU = C, dT and dH = C,,dT, as if these relations 
were applicable to a general change of state, which of 
course they are not. Occasionally the authors themselves 
seem to be led astray in a similar respect. Thus on p. 149- 
150 it is said that the equation 


= C,dT —R 


dS dP 


T P 


is completely general and applicable to all systems. 
Although there has thus been a tendency to over- 

simplify the theory, the book has some very good features 

and most of the writing is clear and straightforward. 


94 


The reviewer enjoyed in particular the chapter on steam- 
power plants. 


The problem of the presentation of thermodynamics 
is one of continual interest. On the one hand, the young 
engineer is usually much more interested in applications 
than in theory and therefore prefers a simplified approach; 
on the other hand, if he fails to achieve a real grasp of the 
fundamentals he is at a loss to know how to proceed 
when he is faced with unfamiliar problems, or when he 
should be making new and inventive applications on his 
own account. Thermodynamics, as we all know, is not 
an easy subject ! K. G. Densicu 


Thermodynamic Functions of Gases. Edited by 
F. Dix, Butterworths Scientific Publications, 1956. 
Vol 1175p. + vi, Vol. IT, 201 p. + vi. 63s. each volume. 


Tue publication of the first two volumes of this new series 
is a very welcome event. Their origin and purpose is 
described by Professor Newrrr, in his preface, in the 
following words “ The preparation of the surveys of the 
thermodynamic properties of gases described in these 
volumes has been undertaken in connection with a scheme 
sponsored by the Department of Scientific and Industrial 
Research in Great Britain. The work is organised by the 
Thermodynamics Committee of the Mechanical Engineer- 
ing Research Board and is being shared among a number 
of academic and industrial organisations which have 
volunteered assistance. The purpose of this scheme is to 
provide reliable thermodynamic data for all industrially 
important gases, and now that the first eight surveys are 
ready, the Committee feels that they should be made 
generally available. Further surveys will follow in due 
course.” 


Each volume contains three or more sections, each 
dealing with a particular gas. Although these various 
sections have been prepared by different authors, they all 
follow a similar pattern, namely (1) survey of existing 
data, (2) construction of the temperature — entropy 
diagram, (3) discussion, (4) tabulation of the thermody- 
namic functions. The procedure agreed between the 
authors was thus to make a critical survey of the existing 
information and then from this data to construct the most 
reliable T — S diagram, including the two-phase region, 
and finally to tabulate values of the entropy, enthalpy, 
specific heat and volume as functions of temperature 
and pressure. 
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The difficult work of editorship has been very ably 
carried out by Dr.* F. Din and he also contributes a 
valuable critical discussion in the first volume on the 
methods of preparing thermodynamic diagrams and tables 
from expeiimental data. As he remarks, the engineer 
would be in a much better position to make the necessary 
allowances when using thermodynamic tables if he realized 
the approximations on which they are based. Tables 
for steam aad other gases are certainly very far from 
revealed truth. 


Thus, in this introductory section, Dr. Din surveys 
the method of constructing the T — S diagram by use 
of P — V — 7 data, together with measurements of the 
heat capacity and the Joule-Thomson coefficient. The 
section is generally very lucid, but a few of the paragraphs 
left the reviewer not quite clear about the author's 
meaning. 


The remaining sections of the first volume, all carefully 
carried through as far as can be seen, deal with ammonia 
(P. Davies), carbon dioxide (D. M. Newrrt, M. U. Pat, 
N. R. Kutoor and J. A. W. Hvueoert) and carbon 
monoxide (SMEETON Lean). The second volume comprises 
similar surveys on air (F. Din), Acetylene (F. Dry), 
ethylene (W. F. L. Dick and A. G. M. Hepiey), propane 
(N. R. Kuvoor and D. M. Newirrr and J. 8S. Bareman) 
and argon (F. Dry). 


Further surveys over another dozen or more gases are 
promised. The volumes will surely become an essential 
addition to any chemical engineering library. 


KENNETH DENBIGH 


F. H. Keatine : Chromium-Nickel Austenitic Steels. 
Butterworths Scientific Publications, London 1956. 
x + 138 pp. 25s. net. 


Tus handbook on stainless steels is the outcome of the 
experience of long usage by the Imperial Chemical 
Industries of these materials. In order to ensure 
that each division of the company is fully informed on 
various subjects these handbooks were initiated for 
internal circulation about twenty years ago. They 
therefore presume text-book knowledge of the subject 
and are designed to bridge the gulf between text-books 
and practice. For stainless steels this is particularly 
valuable since, owing to the large variations in composition, 
it is very difficult for the user to make a fair assessment 
between the claims of various trade materials drawn 
to his attention by various makers. 


The first chapter gives a 6 page introduction to the 
metallurgy of the subject and is followed by 22 pages 
on the historical development of the material. Chapter 
three has a small portion on manufacturing processes 
and a much larger section on various fabricating processes. 
Chapter four is 3 pages on castings. 


The next two chapters on mechanical properties and 
corrosion are most valuable to the chemical engineer 
and discuss the strength of the materials up to 900°C 
and down to — 200°C. A guide is given to suitable 
maximum allowable working stresses for every 100°C 
from room temperature to 1,000°C. The situation regard- 
ing sub-zero temperatures however is still not very clearly 
cut. 


A most useful table of the corrosion resistance of 
18 : 8 steels stabilized with titanium or niobium is given 
for over 100 different chemicals and data is added for the 
effect of additions of molybdenum and copper. 


A chapter on chemical analysis covers carbon, chronium, 
nickel, manganese, silicon, niobium, tungsten, molyb- 
denum, sulphur, titanium and phosphorus. 


The book should go a long way to secure that the next 
review of British Standards on these steels will (to quote 
the author's final words) * provide simple documents 
to cover the purchases of British and foreign users’ and 
not leave important matters of composition to be settled 
at the discretion of the manufacturer. 


M. B. DonaLp 


Physikalisch-chemisches Rechnen in Wissenschaf) 
und Technik. Hans Fromuerz: Verlag Chemie, Wein- 
heim, 1956. xii + 316 pp. DM 32-50 (£2 16s. 11d.) 


Ir is the author’s intention in this book to give instructions 
for the use of physical chemistry in science and technology. 
A large part (206 pages) of some fourteen chapters, is de- 
voted to worked-out problems which, in each chapter, are 
preceded by a summary of the formulae to be used. The 
gas laws, solutions, vapour pressure, thermochemistry, 
chemical equilibrium, electrochemistry, and _ reaction 
kinetics are discussed. Naturally, the author has had to 
confine himself to a selection from the wide range of 
subjects existing in this field. The choice is on the whole 
a good one, but a few criticisms may be made. 


The discussion of the problems of true gases (Chapter 
III), where only the van der Waals equation is dealt with, 
should have been supplemented by generalized state dia- 
grams (particularly in view of their successful application 
in the petroleum industry ). 


In the questions about the determination of molecular 
weights in solutions, etc. (Chapter VII) the activity 
coefficients of the solvent and the solute are not taken into 
consideration, with the result that the determinations only 
apply to strongly diluted solutions. Some problems about 
the determination of activity coefficients would not have 
been out of place in this chapter. 


The excellent Antoine equation could have been men- 
tioned in connection with temperature dependence of 
vapour pressure (Chapter VIII). 





Book reviews 


The chapters on thermochemistry and position of 
equilibrium (Chapters IX, X, XI) are very well handled. 
In these chapters the author should have drawn attention 
of his readers to the existence of the excellent modern 
American tables which give a good summary of a large 
part of the available thermodynamic quantities. 


The last three chapters deal with electrochemistry and 
the reaction kinetics, and in them a good selection of 
problems is given. 


A good symbol and subject index as well as some tables 


with fundamental quantities and conversion factors ter- 
minate the book. 


The foregoing remarks, which partly arise from a per- 
sonal preference, are not meant in any way to disparage 
the high quality of this well-presented book. This work can 
thus be recommended without reservation to chemists and 
chemical engineers, and especially to students of the 
subjects. This reviewer has enjoyed reading the book. 


H, A. G. CaERMIN 





